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ABSTRACT 


The past few years have witnessed a rapid growth in the number and variety of 
applications of fuzzy logic. The applications range from consumer products such as 
cameras, camcorders, washing machines, and microwave ovens to industrial process 
control, medical instrumentation, decision-support systems, and portfolio selection. 

The basic concept underlying fuzzy logic is that of a linguistic variable, that is, a 
variable whose values are words rather than numbers. Although words are inherently less 
precise than numbers, their use is closer to human intuition. Furthermore, computing with 
words exploits the tolerance for imprecision and thereby lowers the cost of solution. 

Another basic concept in fuzzy logic is that of a fuzzy if-then rule or, 
simply, fuzzy rule. The guiding principle of soft computing is to exploit the tolerance for 
imprecision, uncertainty, and partial truth to achieve tractability, robustness, and low cost 
solution. Among various combinations of methodologies in soft computing, the one that 
has highest visibility is that of fuzzy logic. 

Hence the fuzzy logic toolbox of Matlab has been used to develop the fuzzy 
inference system for manufacture of converter lining bricks for Rourkela Steel Plant. 
Data mining tools such as CART (Classification and regression trees) has been used for 
rule generation from data. Rule generation has also been done using ID3. Validation of 
rules by use of Bayesian network and probabilistic study of effect of various parameters 
to predict final output has been done. Two systems have been developed one based on 
data and the other using expert advice for the rule generation. 
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CHAPTER 1 


INTRODUCTION 

Basic Oxygen Furnace (BOF), commonly known as Converter, is a vital link 
in a steel plant and the plant productivity depends largely on the trouble free run of 
the steel-making vessel. Magnesia carbon bricks (MgO-C) are used as lining for the 
converters, that are used for melting of the iron ore and other raw material used in the 
manufacture of steel. Due to the extreme temperatures at which the converters are 
operating (almost 2000 degree centigrade) the lining deteriorates after a specific 
number of manufacturing cycles. 

MgO-C bricks are manufactured at Rourkela Steel Plant in a joint project with 
Research and Development Centre for Iron and Steel (RDCIS), Ranchi. At present 
SAIL engineers have been able to achieve a maximum of up to 1100 cycles of 
converter operation in a single campaign, after which the MgO-C brick lining has to 
be changed leading to a down time of 10 to 15 days. By increasing the life of the 
bricks the number of manufacturing cycles per campaign can be increased, there by 
reducing the frequency of converter lining changes, and hence increasing the 
productivity and reducing cost and saving on man hours. With the globalisation of the 
market the public sector organisation like SAIL have come under excess pressure to 
optimise their production process or be left behind in the competition for market 
share. 

1.1 Problem statement 

To develop a fuzzy inference system model for refractory brick 
manufacturing process at Rourkela steel plant. 

This process involved gathering the vast knowledge and experience of the 
engineers, workers and plant operators. Collection of online data during the 
manufacturing process to substantiate the knowledge acquired and analysing the 
manufacturing process to determine the parameters that have been taken into 
consideration for manufacture of high quality bricks. 

1.2 Organisation of Thesis 

Since a thorough understanding of the under lying process for manufacture of 
the bricks is very important to design the system, chapter 2 has been devoted entirely 
to explain the process in detail. The parameters that critically affect the brick 
properties have been broadly identified. In addition the properties of the raw materials 
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used and the optimised norms to be followed have been tabulated, which gives a 
deeper insight into understanding of the system. 

The software tools used to implement the system are WIN-PROLOG of Logic 
Programming Associates Ltd and Fuzzy Logic toolbox of Matlab. The use of fuzzy 
logic toolbox requires a good foundation of fuzzy logic and fuzzy inference systems 
the salient feature of which shall be discussed in chapter 3. 

Data mining tool that has been used to analyse the data collected from the 
plant for generation of rules are IDS and Cart4, which have been discussed in detail in 
chapter 4 and chapter 5. 

Bayesian network has been used for validation of rules and to carryout a 
probabilistic study of the effect of variation in input parameters to predict final output, 
the details of which have been discussed in chapter 6. 

The fuzzy inference system developed using the vast knowledge acquired 
from the experts working on the plant, right from the engineers to the workers and 
plant operators, is discussed in chapter 7. The fuzzy inference system developed using 
the data collected from the plant is discussed in chapter 8. 

Comparison between the Matlab system and WIN-PROLOG system has been 
brought out in chapter 9. 
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CHAPTER 2 


REFRACTORY BRICK MANUFACTURING PROCESS 

2.1 Basic Concept, brick composition and purpose 

Basic Oxygen Furnace (BOF), commonly known as Converter, is a vital link in a 
steel plant and the plant productivity depends largely on the trouble free run of the steel- 
making vessel. Magnesia carbon bricks (MgO-C) are used as lining for the converters, 
that are used for melting of the iron ore and other raw material used in the manufacture of 
steel. Due to the extreme temperatures at which the converters are operating (almost 2000 
degree centigrade) the lining deteriorates after a specific number of manufacturing cycles. 
At present SAIL engineers have been able to achieve a maximum of up to 1 1 00 cycles of 
converter operation in a single campaign, after which the lining has to be changed leading 
to a down time of 10 to 15 days. By increasing the life of the bricks the number of 
manufacturing cycles can be increased, there by reducing the frequency of converter 
lining changes, and hence increasing the productivity and reducing cost and saving on 
man hours. The purpose of developing the expert system is to optimise the manufacturing 
process so as to steadily achieve higher and higher converter operation life. 

Magnesia carbon bricks essentially constitute of magnesia grains and natural 
graphite. But at room temperature no reaction take place between the two, hence in order 
to manufacture bricks, liquid binder (pitch or resin) is required. At room temperature this 
binder provides the strength. At higher temperature, carbon network is produced from 
this binder and provides necessary binding. 

Manufacture of MgO-C bricks involves preparation of batch with optimum 
granulometry, mixing of these two main components in particular form with liquid 
binder, compaction under appropriate pressure to get optimum bulk density and 
subsequent heat treatment. 

The critical thermo-mechanical properties of MgO-C bricks are largely controlled 
by graphite. Graphite is also responsible for high slag corrosion resistance. Flake 
structure and grain size distribution of the graphite is of decisive importance for the 
press-ability and compaction of the MgO-C bricks as well as oxidation resistance. 
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2.2 The Brick Making Process 

2.2.1 Raw material 

Seawater magnesia in different size fractions, specification of which are 
given in table2.1, along with Natural graphite, specification of which is given in 
table2.2. Pitch is used as liquid binder whose characteristics are given in table2.3. 

2.2.2 Brick making 

Seawater magnesia is first crushed into different grain sizes and filtered 
using screen sieves into the proper granulometry ratio as given in table2.4. This 
weighted quantity of magnesia grains (coarse, medium and fine) are collected in a 
scale car and discharged into the grain heaters, each batch weight is fixed at 1200 
kgs. At start of the operation coarse and medium grains are charged into the grain 
heater drum and heated for 2-4 minutes. The heated grain is then discharged into 
the mixer where the fines along with the graphite are added. Hot mixing is done 
with molten pitch having temperature of about 160 degrees centigrade. The brick 
pressing is done by adjusting the hydraulic pressure to 1450 kg/sqcm. 


Coarse grains Grain heater drum medium grains 



Fig 2.1 brick manufacturing process 



Here the temperature of the mix before pressing is very critical as at 
temperatures exceeding 1 10 degrees centigrade the bricks start cracking, but for 
uniform mixing to be achieved the mixing has to be done at higher temperatures 
of up to 120 to 140 degrees centigrade. Thus we have two contradicting but 
critical requirements to be met which was done by keeping the mix in boxes for 4 
to 8 hours depending on the mix temperature and then were slowly dropped from 
height in to the press hopper and discharged in to another box by rolling back the 
conveyor (which is termed as rolls). Finally the bricks are then sent for tempering 
and then are shrink foil packed. The tree structure of the entire process is shown 
in the fig 2.1. 

2.2.3 Standard Optimised Operating Norms 

The standardised operating norms for the various parameters that have to 
be strictly observed are given in table2.5. 

2.3 Parameters affecting the brick properties. 

The parameters that have been taken into consideration for manufacture of high 
quality bricks are as follows 

• Granulometry 

• Combustion chamber temperature 

• Grain heater drum temperature 

• Grain retention time 

• Mixer temperature 

• Mixing time 

• Mix temperature 

• Storage time 

• Rolls 

• Pitch quantity 

• Graphite quantity 

• Temperature of mix before pressing 

• Forming pressure 

• Green bulk density (GBD) 

The green bulk density thus was the dependent variable, being affected by the 
preceding 13 independent variables, but on closer interaction with the plant engineers and 
the operators it was identified that the granulometry, graphite and pitch were kept 
constant for all campaigns and hence could be ignored. It would only be considered in 
case the output variation was far too high or weird. The GBD was dependent only on 10 
inputs. 
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Figure 2.2: Flowchart of Brick Manufacturing Process TlOl 


variables, but on closer interaction with the plant engineers and the operators it was 
identified that the granulometry, graphite and pitch were kept constant for all campaigns 
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and hence could be ignored. It would only be considered in case the output variation was 
far too high or weird. The GBD was dependent only on 10 inputs. 

The Green Bulk Density is the single most important parameter to be monitored 
because all other parameters in the post manufacture stage are directly or indirectly 
depended on it. Once the MgO_C brick has been manufactured it is sent for tempering, 
after which it is the bricks are warped in polypacks and stored in controlled atmosphere. 
The various factors affecting GBD and the parameters that are measured to check the 
quality of the bricks is enumerated below, 

In the plant the important parameters which affect the Green Bulk Density are 
listed below: 

• Combustion Chamber Temperature 

• Grain Heater Drum Temperature 

• Retention Time in Grain Heater 

• Grain Temperature 

• Mixer temperature 

• Mixing time 

• Mix temperature 

• Storage time 

• Number of rolls 

• Press Mix Temperature 

• Forming pressure 

• Pitch quantity 

• Graphite quantity 

The important properties of the MgO-C bricks before and after tempering are 

listed below: 

• Green Bulk Density 

• Tempered Bulk Density 

• Cold crushing strength 

• Coked porosity 

• Hot Modulus of rupture 

• Slag corrosion index 
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2.3.1 Grain Temp (GT) 

Grain temperature is the temperature of the coarse grains when being fed 
to the mixer. In the plant their does not exists any direct means to measure this 
temperature. The temperature has to be measured manually by using a thermal 
sensor attached to a rod by pushing it deep into the heap of grains. This 
temperature is given importance due to the very fact that it affects the temperature 
of mix and its quality directly. Also the heat transfer from coarse grains to the fine 
grains and powdered graphite takes place. The grain temperature depends on the 
following parameters: 

♦♦♦ Combustion Chamber Temperature Initial (CTI) 

❖ Grain Heater Drum Temperature (GHDT) 

♦♦♦ Retention Time in Grain Heater (RT) 

The more is the value of combustion chamber temperature, grain heater drum or 
retention time; the higher will be the grain temperature. 

2.3.2 Mix Temp (MT) 

Mix Temperature is the temperature of the fully prepared mixture after 
mixing of all the grain types, graphite and dosing of liquid pitch. If the 
temperature is obtained as per specification then the mix can directly taken for 
making of the bricks in the plant press. The Mix temperature depends on 
following parameters: 

❖ Grain Temperature (GT) 

❖ Mixer temperature (Mix_Temp) 

❖ Mixing time 

Higher the grain temperature, Mixer temperature and the Mixing Time, 
higher will be the mix temperature. 

2.3.3 Press Mix Temperature (PMT) 

Press Mix Temperature is the Mix Temperature at the time of the mix 
going in for making of the bricks. It happens very regularly that the Mix 
temperature varies from its actual required value and the mix has to be cooled 
down or reheated by repeating the process of heating. Cooling takes place by just 
leaving the drum containing mix on one side for some time and let it cool 
naturally. This is called Storage time. Other method to cool the mix quickly is to 
raise the drum to certain height and pour the mix slowly into another drum on 
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ground thus letting it dissipate heat to surrounding air. One such process is called 
rolling of the mix and this process can be repeated a number of times. Thus the 
Press Mix Temperature will depend on following parameters: 

> Mix temperature 

> Storage time 

> Number of rolls 

The Mix temperature is brought to the specified value of Press Mix 
temperature by storage time and the number of rolls. All temperatures in the plant 
are measured in degrees centigrade and hence for obvious reasons the units have 
not been mentioned explicitly any where. 

2.3.4 Green Bulk Density (GBD) 

This is the most important property of brick in the plant process. If the 
green BD is achieved correctly, the other properties given the correct composition 
of materials and their quality will generally will be as per the specifications. It is 
measured in grams per cubic centimeters and generally varies from 2.92 to 3.15. 
Generally the higher the density of the brick greater is the life achieved. Green 
BD will directly depend on the following parameters: 

> Press Mix Temperature 

> Forming pressure 

> Pitch quantity 

> Graphite quantity 

2.3.5 Cold Crushing Strength (CCS) 

It is the measure of the strength required to rupture the brick under 
compression at room temperature. It is measured in kilograms per square 
centimeter. Higher the GBD and Pitch higher will be the CCS, but the Tempering 
should to optimized, its value varies between 212 to 354 Kg/cm 2 . The CCS of 
bricks depends on following: 

• Pitch quantity 

• Green Bulk Density 

• Existing Tempering 



2.3.6 Coked Porosity (CP) 

It varies depending on the quantity of retained carbon and its value ranges 
from 10.82 to 15.25 %, thus it is measured as a % of the retained carbon. The CP 
of bricks depends on the following: 

■ Pitch quantity 

■ Green Bulk Density 

2.3.7 Hot Modulus of rupture (HMOR) 

Modulus of rupture is the maximum stress that a rectangular test piece of 
specified dimension can withstand when it is bent in a three point bending device. 
It is tested in both oxidizing atmosphere and reducing atmosphere, it a measure of 
the strength of the bricks under heated condition. The test pieces are heated to the 
test temperature (generally 1400°C for MgO-C). It has a higher value under 
reducing condition that oxidizing condition. The typical values under reducing 
condition vary between 60 and 105 Kg/cm 2 , under oxidizing condition they vary 
between 15 and 35. The HMOR depends on the following: 

❖ Green Bulk Density 

❖ Cold crushing strength 

❖ Carbon Content 

❖ Aluminum powder 

2.3.8 Slag Corrosion Resistance (SCR) 

It generally depends on the quantity of graphite added, bricks with higher 
graphite showed lower SCR and those with lower graphite showed higher SCR. 
This is also specified in terms of Slag Corrosion index and measured as the % of 
wear in millimeters. Its typical values vary between 100 and 120 %. The SCR 
depends on the following parameters: 

❖ Carbon Content 

❖ Chemical purity 

❖ Coked Porosity 
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Table2.1 Specification of sea water magnesia used 


Properties chemical 

composition % 

Source 

Sardamag Italy 

Periclase, Israel 

Periclase, India 

MgO min 

97 

99 

97 

CaO max 

2 

0.8 

2 

SiC >2 max 

0.5 

0.05 

0.35 

FeCb max 

0.5 

0.07 

0.18 

AI 2 O 3 max 

— 

0.05 

0.12 

B 2 O 3 max 

0.15 

0.005 

0.03 

Bulk density (gm/cc) 

— 

3.42 

3.4 

Crystal size {pi) 

— 

70-90 

100 


Table2.2 Typical Granulometry ratio 


Grain size( mm) 

Norm (%) 

5-8 

15 

3-5 

20 

1.6-3 

20 

0.5-1 . 6 

15 

0-0.5 

10 

0 - 0.1 

20 


Table2.3 Typical Optimised Operating Parameter Norms 


Parameters 

Norms 

Granulometry 

Tolerance +/- 10 kg 

Pitch quantity 

60 +/- 5 kg 

Graphite quantity 

60 +/- 2 kg 

Mix temperature 

120+/- 10 U C 

Pressing temperature 

105+/- 5 °C 

De-airing 

Minimum twice 

Green bulk density 

> 3 .00 gm/cc 

CCS (tempered) 

> 350 kg/cm 2 

Bulk density (tempered) 

> 2.98 gm/cc 
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CHAPTER 3 


FOUNDATIONS OF FUZZY LOGIC AND FUZZY INFERENCE 

SYSTEMS 

3.1 Basic Concepts: - Fuzzy Logic 

Fuzzy logic is all about the relative importance of precision: How important is it 
to be exactly right when a rough answer will do? Fuzzy logic is best summarised no 
better than in the words of the two great scientists, 

So far as the laws of mathematics refer to reality, they are not certain. And so far 
as they are certain, they do not refer to reality. — Albert Einstein 

As complexity rises, precise statements lose meaning and meaningful statements 
lose precision. — Lotfi Zadeh 

Fuzzy logic is a convenient way to map an input space to an output space. This is 
the starting point for everything else, and the great emphasis here is on the word 
“convenient.” What do you mean by mapping input space to output space? Here are a few 
examples: You tell me how good your service was at a restaurant, and I will tell 
You what the tip should be. You tell me how hot you want the water, and I will adjust the 
faucet valve to the right setting. You tell me how far away the subject of your photograph 
is, and 1 will focus the lens for you. You tell me how fast the car is going and how hard 
the motor is working, and I will shift the gears for you. 


Input Space Output Space 

(all possible service (ail possible tips) 

quality ratings) 



Fie 3.1 An input-output map for the tipping problem: 

It’s all just a matter of mapping inputs to the appropriate outputs (Fig 3.1). Between the 
input and the output we’ll put a black box that does the work. What could go in the black 
box? Any number of things: fuzzy systems, linear systems, expert systems, neural 
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networks, differential equations, interpolated multidimensional lookup tables, or even a 
spiritual advisor, just to name a few of the possible options. Clearly the list could go on 
and on. 

3.2 Why Fuzzy and not others 

Of the dozens of ways to make the black box shown above work, it turns out that 
fuzzy is often the very best way. Why should that be? As Lotfi Zadeh, who is considered 
to be the father of fuzzy logic, once remarked: “In almost every case you can build the 
same product without fuzzy logic, but fuzzy is faster and cheaper.” 

Listed below are some of the advantages offered by fuzzy logic 
•Fuzzy logic is conceptually easy to understand. 

The mathematical concepts behind fuzzy reasoning are very simple. What makes 
fuzzy nice is the “naturalness” of its approach and not its far-reaching complexity. 
•Fuzzy logic is flexible. 

With any given system, it’s easy to massage it or layer more functionality on top 
of it without starting again from scratch. 

•Fuzzy logic is tolerant of imprecise data. 

Everything is imprecise if you look closely enough, but more than that, most 
things are imprecise even on careful inspection. Fuzzy reasoning builds this 
understanding into the process rather than tacking it onto the end. 

•Fuzzy logic can model nonlinear functions of arbitrary complexity. 

You can create a fuzzy system to match any set of input-output data. This process 
is made particularly easy by adaptive techniques like ANFIS (Adaptive Neuro- 
Fuzzy Inference Systems), which are available in the Fuzzy Logic Toolbox. 

•Fuzzy logic can be built on top of the experience of experts. 

In direct contrast to neural networks, which take training data and generate 
opaque, impenetrable models, fuzzy logic lets you rely on the experience of 
people who already understand your system. 

•Fuzzy logic can be blended with conventional control techniques. 

Fuzzy systems don’t necessarily replace conventional control methods. In many 
cases fuzzy systems augment them and simplify their implementation. 

•Fuzzy logic is based on natural language. 

The basis for fuzz y logic is the basis for human communication. This observation 
underpins many of the other statements about fuzzy logic. 
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The last statement is perhaps the most important one and deserves more discussion. 
Natural language, that which is used by ordinary people on a daily basis, has been shaped 
by thousands of years of human history to be convenient and efficient. Sentences written 
in ordinary language represent a triumph of efficient communication. We are generally 
unaware of this because ordinary language is, of course, something we use every day. 
Since fuzzy logic is built atop the structures of qualitative description used in everyday 
language, fuzzy logic is easy to use. 

3.3 When Not to Use Fuzzy Logic 

Fuzzy logic is not a cure-all. When should you not use fuzzy logic? The safest 
statement is the first one made in this introduction: fuzzy logic is a convenient way to 
map an input space to an output space. If you find it’s not convenient, try something else. 
If a simpler solution already exists, use it. Fuzzy logic is the codification of common 
sense; use common sense when you implement it and you will probably make the right 
decision. Many controllers, for example, do a fine job without using fuzzy logic. 
However, if you take the time to become familiar with fuzzy logic, you’ll see it can be a 
very powerful tool for dealing quickly and efficiently with imprecision and nonlinearity. 

3.4 Foundations of Fuzzy Logic 

Everything is vague to a degree you do not realize till you have tried to make it 
precise. - Bertrand Russell 

3.4.1 Fuzzy Sets 

Fuzzy logic starts with the concept of a fuzzy set. A fuzzy set is a set 
without a crisp, clearly defined boundary. It can contain elements with only a 
partial degree of membership. To understand what a fuzzy set is, first consider 
what is meant by what we might call a classical set. A classical set is a container 
that wholly includes or wholly excludes any given element. For example, the set 
of days of the week unquestionably includes Monday, Thursday, and Saturday. It 
just as unquestionably excludes butter, liberty, and dorsal fins, and so on. We call 
this set a classical set because it’s been around for such a long time. It was 
Aristotle who first formulated the Law of the Excluded Middle, which says X 
must either be in set A or in set not A. Another version runs like this. Of any 
subject, one thing must be either asserted or denied. 

Here is a restatement of the law with annotations: “Of any subject (say 
Monday), one thing (being a day of the week) must be either asserted or denied (I 
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assert that Monday is a day of the week).” This law demands that opposites, the 
two categories A and not-A, should between them contain the entire universe. 
Everything falls into either one group or the other. There is no thing that is both a 
day of the week and not a day of the week. 

Now consider the set of days comprising a weekend. The fig 3.2 below is 
one attempt at classifying the weekend days. 



Fig 3.2 Non-classical set 

Most would agree that Saturday and Sunday belong, but what about 
Friday? It “feels” like a part of the weekend, but somehow it seems like it should 
be technically excluded. So in the diagram above Friday tries its best to sit on the 
fence. Classical or “normal” sets wouldn’t tolerate this kind of thing. Either 
you’re in or you’re out. Human experience suggests something different, though: 
fence sitting is a part of life. 

Of course we’re on tricky ground here, because we’re starting to take 
individual perceptions and cultural background into account when we define what 
constitutes the weekend. But this is exactly the point. Even the dictionary is 
imprecise, defining the weekend as “the period from Friday night or Saturday to 
Monday morning.” We’re entering the realm where sharp edged yes-no logic 
stops being helpful. Fuzzy reasoning becomes valuable exactly when we’re 
talking about how people really perceive the concept “weekend” as opposed to a 
simple-minded classification useful for accounting purposes only. More than 
anything else, the following statement lays the foundations for fuzzy logic. 

In fuzzy logic, the truth of any statement becomes a matter of degree. 

Any statement can be fuzzy. The tool that fuzzy reasoning gives is the 
ability to reply to a yes-no question with a not-quite-yes-or-no answer. This is the 
kind of thing that humans do all the time (think how rarely you get a straight 
answer to a seemingly simple question) but it’s a rather new trick for computers. 
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How does it work? Reasoning in fuzzy logic is just a matter of 
generalizing the familiar yes-no (Boolean) logic. If we give “true” the numerical 
value of 1 and “false” the numerical value of 0, we’re saying that fuzzy logic also 
permits in-between values like 0.2 and 0.7453. For instance, 

Q: Is Saturday a weekend day? 

A: 1 (yes, or true) 

Q: Is Tuesday a weekend day? 

A: 0 (no, or false) 

Q: Is Friday a weekend day? 

A: 0.8 (for the most part yes, but not completely) 

Q: Is Sunday a weekend day? 

A: 0.95 (yes, but not quite as much as Saturday). 

Below in fig 3.5 is a plot that shows the truth-values for “weekend-ness” if 
we are forced to respond with an absolute yes or no response. Fig 3.3 below is a 
plot that shows the truth-value for weekend-ness if we are allowed to respond 
with fuzzy in-between values. 



THUSDAY FRIDAY SATURDAY SUNDAY MONDAY 


FIG 3.3 DAYS OF WEEKED TWO VALUED FUNCTION 

Technically, the representation is from the domain of multivalued logic (or 
multivalent logic). If I ask the question “Is X a member of set A?” the answer might be 
yes, no, or any one of a thousand intermediate values in between as in fig 3.4. In other 
words, X might have partial membership in A. Multivalued logic stands in direct contrast 
to the more familiar concept of two-valued (or bivalent yes-no) logic. To return to our 
example, now consider a continuous scale time plot of weekend-ness shown in fig 3.6. 
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FIG 3.4 DAYS OF WEEKED MULTI-VALUED FUNCTION 


1.0 

Weekend-ness 

0.0 


Thus Fri Sat Sun Mon 

Fig 3.5 Days Of The Weekend Two - Valued Membership 



Fig 3.6 Days Of The Weekend Multi- Valued Membership 

By making the plot continuous, we’re defining the degree to which any given 
instant belongs in the weekend rather than an entire day. In the plot on the left, notice that 
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at midnight on Friday, just as the second hand sweeps past 12, the weekend-ness truth- 
value jumps discontinuously from 0 to 1. This is one way to define the weekend, and 
while it may be useful to an accountant, it doesn’t really connect with our real-world 
experience of weekend-ness. 

One plot shows a smoothly varying curve that accounts for the fact that all of 
Friday, and, to a small degree, parts of Thursday, partake of the quality of weekend-ness 
and thus deserve partial membership in the fuzzy set of weekend moments. The curve 
that defines the weekend-ness of any instant in time is a function that maps the input 
space (time of the week) to the output space (weekend-ness). Specifically it is known as a 
membership function. 

3.4.2 Membership Functions 

A membership function (MF) is a curve that defines how each point in the input 
space is mapped to a membership value (or degree of membership) between 0 and 1 . The 
input space is sometimes referred to as the universe of discourse. 

One of the most commonly used examples of a fuzzy set is the set of tall people. 
In this case the universe of discourse is all potential heights, say from 3 feet to 
9 feet, and the word “tall” would correspond to a curve that defines the degree to which 
any person is tall. If the set of tall people is given the well-defined (crisp) boundary of a 
classical set, we might say all people taller than six feet are officially considered tall. But 
such a distinction is clearly absurd. It may make sense to consider the set of all real 
numbers greater than six because numbers belong on an abstract plane, but when we want 
to talk about real people, it is unreasonable to call one person short and another one tall 
when they differ in height by the width of a hair. 

excellent! 



Fig 3.7 Distinction between Tall and Short 
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1.0 


tall (u = 1.0) 


degree of 
membership, p 


0.0 




sharp-edged 
membership 
function for 
TALL 



not tall (ji = 0.0) 


height 

Fig 3.8 Two valued membership function 

But if the kind of distinction shown in Fig 3.7 and Fig 3.8 above is unworkable, 
then what is the right way to define the set of tall people? Much as with our plot of 
weekend days, the Fig 3.9 below shows a smoothly varying curve that passes from not 
tall to tall. The output-axis is a number known as the membership value between 0 and 1. 
The curve is known as a membership function and is often given the designation of ju . 
This curve defines the transition from not tall to tall. Both people are tall to some degree, 
but one is significantly less tall than the other. 


definitely a tall 
person (u = 0.95) 


really not very 
tall at all (p = 0.30) 


Fig 3.9 Days Of The Weekend Multi- Valued Membership 

Subjective interpretations and appropriate units are built right into fuzzy sets. If I 
say, “She’s tall,” the membership function “tall” should already take into account 
whether I’m referring to a six-year-old or a grown woman. Similarly, the units are 
included in the curve. Certainly it makes no sense to say “Is she tall in inches or in 
meters?” 

3.4.3 Logical Operations 

We now know what’s fuzzy about fuzzy logic, but what about the logic? The 
most important thing to realize about fuzzy logical reasoning is the fact that it is a 
superset of standard Boolean logic. In other words, if we keep the fuzzy values at their 
extremes of 1 (completely true), and 0 (completely false), standard logical operations will 
hold. As an example, consider the standard truth tables in Fig 3.10 below. 


1.0 

conti nu ous 


memberships^ 

degree of 

function for jf 

membership, p 


0.0 

] ' 


height 
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extremes of 1 (completely true), and 0 (completely false), standard logical operations will 
hold. As an example, consider the standard truth tables in Table 3.1 below. 


AND 

OR 

NOT 

A 

B 

A and B 

A 

B 

A or B 

A 

not A 

0 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

1 

1 

1 

0 

1 

0 

0 

1 

0 

1 


1 

1 

1 

1 

1 

1 


Table 3.1 Standard Truth Tables 

Now remembering that in fuzzy logic the truth of any statement is a matter of 
degree, how will these truth tables be altered? The input values can be real numbers 
between 0 and 1. What function will preserve the results of the AND truth table (for 
example) and also extend to all real numbers between 0 and 1? One answer is the min 
operation. That is, resolve the statement A AND B, where A and B are limited to the range 
(0,1), by using the function min(A,B). Using the same reasoning, we can replace the OR 
operation with the max function, so that A OR B becomes equivalent to max (A, B). 
Finally, the operation NOT A becomes equivalent to the operation. Notice how the truth 
table above is completely unchanged by this substitution in Table 3.2 below. 


Min 

Max 

Compliment 

A 

B 

Min(A,B) 

A 

B 

Max(A,B) 

A 

1-A 

0 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

1 

1 

1 

0 

1 

0 

0 

1 

0 

1 


1 

1 

1 

1 

1 

1 


Table 3.2 Standard Truth Tables equivalent to min max and compliment 


Moreover, since there is a function behind the truth table rather than just the truth 
table itself, we can now consider values other than 1 and O.Fig 3.12 uses a graph to show 
the same information. We have converted the truth table to a plot of two fuzzy sets 
applied together to create one fuzzy set. The upper part of the figure displays plots 
corresponding to the two-valued truth tables above, while the lower part of the figure 
displays how the operations work over a continuously varying range of truth values A and 
B according to the fuzzy operations we have defined. 
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AND OR NOT 

maxfA.B) fl-A] 


Fig 3.10 Graphical representation of Standard Operators 

Given these three functions, we can resolve any construction using fuzzy sets and 
the fuzzy logical operation AND, OR, and NOT. 

3.4.4 If-Then Rules 

Fuzzy sets and fuzzy operators are the subjects and verbs of fuzzy logic [9], These 
if-then rule statements are used to formulate the conditional statements that comprise 
fuzzy logic. 

A single fuzzy if-then rule assumes the form 
ifx isX theny is B 

Where A and B are linguistic values defined by fuzzy sets on the ranges (universes 
of discourse) X and Y, respectively. The if-part of the rule “x is X” is called the 
antecedent or premise, while the then-part of the rule “y is 5” is called the consequent or 
conclusion. An example of such a rule might be 

if service is good then tip is average 

Note that good is represented as a number between 0 and 1, and so the antecedent 
is an interpretation that returns a single number between 0 and 1. On the other hand, 
average is represented as a fuzzy set, and so the consequent is an assignment that assigns 
the entire fuzzy set B to the output variable y. In the if-then rule, the word “is” gets used 
in two entirely different ways depending on whether it appears in the antecedent or the 
consequent. A less confusing way of writing the rule would be 

if service == good then tip = average 
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In general, the input to an if-then rule is the current value forthe input variabDeQn 
this case, service ) and the output is an entire fuzzy set (in this case, averages). "Thilss=et 
will later be defuzzified, assigning one value to the output. 

Interpreting an if-then rule involves distinct parts: first evaluating t he antecede nt 
(which involves fuzzifying the input and applying any necessary fiezzy Ojgxr-alcm) • and 
second applying that result to the consequent (known as implication). In the case of tw«c- 
valued or binary logic, if-then rules don’t present much difficulty. If the premise is trine, 
then the conclusion is true. If we relax the restrictions of two-valued log ic and le _ t the 
antecedent be a fuzzy statement, how does this reflect on the conclusion? Tbieanswrewis a 
simple one. If the antecedent is true to some degree of membership, then thie com seqpemt 
is also true to that same degree. In other words 

in binary logic: p -> q(p and q are either both true or both false) 

in fuzzy logic: 0.5p->0.5 q (partial antecedents provide partial implicati on)) 

The antecedent of a rule can have multiple parts. 

if sky is gray and -wind is strong and barometer is falling, then . .. 
in which case all parts of the antecedent are calculated simultaneously and resolwed to a 
single number using the logical operators described in the preceding section, Thie 
consequent of a rule can also have multiple parts. 

if temperature is cold then hot water valve is open and cold- wee ter verhe is ~shwt 
in which case all consequents are affected equally by the result of the antecedent. Blow 
does the antecedent affect the consequent? The consequent specifies a fuzzy se=ttoe 
assigned to the output. The implication function then modifies that fuzzy set to the degree 
specified by the antecedent. The most common ways to modify the output fuzzy sen am 
truncation using the min function (where the fuzzy set is “chopped off’ as shown be low) 
or scaling using the prod function (where the output fuzzy set is c ‘squashe <f>. 

3. 4.4.1 Summary of If-Then Rules 

Interpreting if-then rules is a three-part process. Fig 3. 13 below illustrates these pa. rts. 

(1) Fuzzify inputs : Resolve all fuzzy statements in the antecedent to a degreeecof 
membership between 0 and 1. If there is only one part to the antecedent, thdsissthae 
degree of support for the rule. 

(2) Apply fuzzy operator to multiple part antecedents. If there are multiple parlstto 
the antecedent, apply fuzzy logic operators and resolve the antecedent to a si ngle 
number between 0 and 1. This is the degree of support for the rule . 
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(3) Apply Mmplication method : Use the degree of support for the entire rule to 
shape the output fuzzy set. The consequent of a fuzzy rule assigns an entire fuzzy 
set to the output. This fuzzy set is represented by a membership function that is 
chosen to indicate the qualities of the consequent. If the antecedent is only 
partially true, (i.e., is assigned a value less than 1), then the output fuzzy set is 
truncated according to the implication method. 

In general, one rule by itself doesn't do much good. What are needed are two or more 
rules that can play off one another. The output of each rule is a fuzzy set. The output 
fuzzy sets for each rule are then aggregated into a single output fuzzy set. Finally the 
resulting set is defuzzified , or resolved to a single number. The next section shows how 
the whole process works from beginning to end for a particular type of fuzzy inference 
system called a Mamdani type. 


Antecedent 


Consequent 


1. Fuzzify 
inputs 


2. Apply 
OR operator 
(max) 


il Apply 
implication 
operator (min) 


If service is excellent or food is delicious then tip = generous 



P (sen/ice=~exceiien i) ~ 0 .0 P{food==ddichus) - 0 .7 


( 0.0 


or 


0.7 ) 


then tip = generous 



max(OMpOJ) = OJ 


( 0.7 ) 


then tip - generous 



1 

A 



Q.7 

r 

generous § A 




mm(0J 3 generous) tip (fuzzy) 


Fig 3.13 Interpreting If-Then Rules 3 part process 
3.5 Fuzzy Inference Systems 

Fuzzy inference is the process of formulating the mapping from a given input to 
an output using fuzzy logic. The mapping then provides a basis from which decisions can 
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be made, or patterns discerned. The process of fuzzy inference involves all of the pieces 
that are described in the previous sections: membership functions, fuzzy logic operators, 
and if-then rules. There are two types of fuzzy inference systems that can be implemented 
in the Fuzzy Logic Toolbox: Mamdani-type and Sugeno-type. These two types of 
inference systems vary somewhat in the way outputs are determined. 

Fuzzy inference systems have been successfully applied in fields such as 
automatic control, data classification, decision analysis, expert systems, and computer 
vision. Because of its multidisciplinary nature, fuzzy inference systems are associated 
with a number of names, such as fuzzy-rule-based systems, fuzzy expert systems, fuzzy 
modeling, fuzzy associative memory, fuzzy logic controllers, and simply (and 
ambiguously) fuzzy systems. Since the terms used to describe the various parts of the 
fuzzy inference process are far from standard, we will try to be as clear as possible about 
the different terms introduced in this section. 

Mamdani’s fuzzy inference method is the most commonly seen fuzzy 
methodology. Mamdani’s method was among the first control systems built using fuzzy 
set theory. It was proposed in 1 975 by Ebrahim Mamdani as an attempt to control a steam 
engine and boiler combination by synthesizing a set of linguistic control rules obtained 
from experienced human operators. Mamdani’s effort was based on Lotfi Zadeh’s 1973 
paper on fuzzy algorithms for complex systems and decision processes. Although the 
inference process we describe in the next few sections differs somewhat from the 
methods described in the original paper, the basic idea is much the same. Mamdani-type 
inference expects the output membership functions to be fuzzy sets. After the aggregation 
process, there is a fuzzy set for each output variable that needs defuzzification. 

It’s possible, and in many cases much more efficient, to use a single spike as the 
output membership function rather than a distributed fuzzy set. This is sometimes known 
as a singleton output membership function, and it can be thought of as a pre-defuzzified 
fuzzy set. It enhances the efficiency of the defuzzification process because it greatly 
simplifies the computation required by the more general Mamdani method, which finds 
the centroid of a two-dimensional function. Rather than integrating across the two- 
dimensional function to find the centroid, we use the weighted average of a few data 
points. Sugeno-type systems support this type of model. In general, Sugeno-type systems 
can be used to model any inference system in which the output membership functions are 
either linear or constant. 
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3.5.1 Five Steps Leading to Development of Fuzzy Inference System 
Step 1. Fuzzily Inputs 

The first step is to take the inputs and determine the degree to which they 
belong to each of the appropriate fuzzy sets via membership functions. In the 
Fuzzy Logic Toolbox, the input is always a crisp numerical value limited to the 
universe of discourse of the input variable (in this case the interval between 0 and 
10) and the output is a fuzzy degree of membership in the qualifying linguistic set 
(always the interval between 0 and 1). Fuzzification of the input amounts to either 
a table lookup or a function evaluation. 

The example we are using in this section is built on three rules, and each 
of the rules depends on resolving the inputs into a number of different fuzzy 
linguistic sets: service is poor, service is good, food is rancid, food is delicious, 
and so on. Before the rules can be evaluated, the inputs must be fuzzified 
according to each of these linguistic sets. For example, to what extent is the food 
really delicious? The Fig 3.14 below shows how well the food at our hypothetical 
restaurant (rated on a scale of 0 to 10) qualifies, (via its membership function), as 
the linguistic variable “delicious.” In this case, we rated the food as an 8, which, 
given our graphical definition of delicious, corresponds to p. = 0.7 for the 
“delicious” membership function. 


1. Fuzzily 
inputs. 


food is delicious 

food = 8 

input 

Fig 3.14 Membership Function For Variable Delicious 

(The compliment to the chef would be “your food is delicious to the degree O.7.”) 
In this manner, each input is fuzzified over all the qualifying membership 
functions required by the rules. 



0.7 


Result of 
fuzzification 
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Step 2. Apply Fuzzy Operator 

Once the inputs have been fuzzified, we know the degree to which each 
part of the antecedent has been satisfied for each rule. If the antecedent of a given 
rule has more than one part, the fuzzy operator is applied to obtain one number 
that represents the result of the antecedent for that rule. This number will then be 
applied to the output function. The input to the fuzzy operator is two or more 
membership values from fuzzified input variables. The output is a single truth- 
value. 

As is described in the section on fuzzy logical operations, any number of 
well-defined methods can fill in for the AND operation or the OR operation. In 
the Fuzzy Logic Toolbox, two built-in AND methods are supported: min 
(minimum) and prod (product). Two built-in OR methods are also supported: max 
(maximum), and the probabilistic OR method probor. The probabilistic OR 
method (also known as the algebraic sum) is calculated according to the equation 
probor (a, b) = a + b - ab 

In addition to these built-in methods, you can create your own methods for 
AND and OR by writing any function and setting that to be your method of 
choice. There will be more information on how to do this later. 

Fig 3.15 below is an example of the OR operator max at work. We’re 
evaluating the antecedent of the rule 3 for the tipping calculation. The two 
different pieces of the antecedent (service is excellent and food is delicious) 
yielded the fuzzy membership values 0.0 and 0.7 respectively. The fuzzy OR 
operator simply selects the maximum of the two values, 0.7, and the fuzzy 
operation for rule 3 is complete. If we were using the probabilistic OR method, 
the result would still be 0.7 in this case. 

Step 3. Apply Implication Method 

Before applying the implication method, we must take care of the rule’s 
weight. Every rule has a weight (a number between 0 and 1), which is applied to 
the number given by the antecedent. Generally this weight is 1 (as it is for this 
example) and so it has no effect at all on the implication process. From time to 
time you may want to weight one rule relative to the others by changing its weight 
value to something other than 1 . Once proper weighting has been assigned to each 
rule, the implication method is implemented as in Fig 3.16. 
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1. Fuzzily 
■inputs. 


2, Apply 

OR operator (max}. 



0 3 

result of 
fuzzy operator 


service = 3 


food = 8 


input 2 


Fig 3.15 Probabilistic OR operate r 
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Anteoed «snt 
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Consequent 


■> 


1. Fuzzi fy 3. Apply 

mputs. Off operator (fmx).. Impffcaticm 

operator (mm). 





If service is excellent or ' foo-d Is dell clou, s tlion tip- generous 


result of 
implication 


servlc e = 3 food = B 

input 1 input 2 


Fig 3.16 Im plication Method 

A consequent is a fuzzy set represented by a membership function, which 
weights appropriately the lingui Stic characteristics that are attributed to it. The 
consequent is reshaped using a function associated with the antecedent (a single 
number). The input for the impl ica_tion process is a single number given by the 
antecedent, and the output is a fuzzy set I cmpli cation is implemented for each rule. 
Two built-in methods are supported, andthey are the same functions that are used 
by the AND method: min (m inimam), viiich truncates the output fuzzy set, and 
prod (product), which scales the output fuzzy set. 

Step 4. Aggregate All Outputs 

Since decisions are based on the testing of all of the rules in an FIS, the 
rules must be combined in some manner im order to make a decision. Aggregation 
is the process by which the fuzzy sets that represent the outputs of each rule are 
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combined into a single fuzzy set. Aggregation only occurs once for each output 
variable, just prior to the fifth and final step, defuzzification. The input of the 
aggregation process is the list of truncated output functions returned by the 
implication process for each rule. The output of the aggregation process is one 
fuzzy set for each output variable. 

Notice that as long as the aggregation method is commutative (which it 
always should be), then the order in which the rules are executed is unimportant. 
Three built-in methods are supported: max (maximum), probor (probabilistic or), 
and sum (simply the sum of each rule’s output set). 

In the Fig 3.17 below, all three rules have been placed together to show 
how the output of each rule is combined, or aggregated, into a single fuzzy set 
whose membership function assigns a weighting for every output (tip) value. 


2. Appfy 3. Appty 

?. Fuzzify inputs. fuzzy rmpmmm 

coBfsHort msthocf frnnrrj. 



If service is poor or food is rancid then tip = cheap 





service « 3 food « 8 

Input 1 input 2 



aggregation 


Fig 3.17 Aggregation Method 
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Step 5. Defuzzify 

The input for the defuzzification process is a fuzzy set (the aggregate 
output fuzzy set) and the output is a single number. As much as fuzziness helps 
the rule evaluation during the intermediate steps, the final desired output for each 
variable is generally a single number. However, the aggregate of a fuzzy set 
encompasses a range of output values, and so must be defuzzified in order to 
resolve a single output value from the set. 

Perhaps the most popular defuzzification method is the centroid 
calculation, which returns the center of area under the curve as shown in Fig 3.18 
below. There are five built-in methods supported: centroid, bisector, middle of 
maximum (the average of the maximum value of the output set), largest of 
maximum, and smallest of maximum. 


5. Defuzzify the 
aggregate output 
(centroid)* 


Fig 3.18 Centriod Method for Defuzzification 
3.5.2 The Fuzzy Inference Diagram 

The fuzzy inference diagram in Fig 3.19 below is the composite of all the smaller 
diagrams we have been looking at so far in this section. It simultaneously displays all 
parts of the fuzzy inference process we have examined. Information flows through the 
fuzzy inference diagram as shown below. 

Notice how the flow proceeds up from the inputs in the lower left, then across 
each row, or rule, and then down the rule outputs to finish in the lower right. This is a 
very compact way of showing everything at once, from linguistic variable fuzzification 
all the way through defuzzification of the aggregate output. 



tip = 16.7% 

Result of 
defuzzification 
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Interpreting the 
Fuzzy Inference 
Diagram 



Fig 3.19 Fuzzy Inference Diagram 
3.5.3Sugeno-Type Fuzzy Iuference 

The fuzzy inference process we’ve been referring to so far is known as 
Mamdani’s fuzzy inference method. It’s the most commonly seen fuzzy methodology. In 
this section we discuss the so-called Sugeno, method of fuzzy inference. It is similar to 
the Mamdani method in many respects. In fact the first two parts of the fuzzy inference 
process, fuzzifying the inputs and applying the fuzzy operator, are exactly the same. The 
main difference between Mamdani type of fuzzy inference and Sugeno-type is that the 
output membership functions are only linear or constant for Sugeno-type fuzzy inference. 

A typical fuzzy rule in a zero-order Sugeno fuzzy model has the form 

if x is A and y is B then z = k 

Where A and B are fuzzy sets in the antecedent, while k is a crisply defined 
constant in the consequent. When the output of each rule is a constant like this, the 
similarity with Mamdani’s method is striking. The only distinctions are the fact that all 
output membership functions are singleton spikes, and the implication and aggregation 
methods are fixed and cannot be edited. The implication method is simply multiplication, 
and the aggregation operator just includes all of the singletons. 

The Fig 3.20 below shows the fuzzy tipping model developed in previous sections 
of this manual adapted for use as a zero-order Sugeno system. Fortunately it is frequently 
the case that singleton output functions are completely sufficient for a given problem’s 
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Fig 3.20 Sugeno fuzzy system for the tipping problem 

The Mamdani implication method has been used in design of the Expert based 
system for RSP and the Sugeno method has been used in design of the data based system. 
Hence this has been discussed in detail in this chapter. 
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CHAPTER 4 


DECISION TREES (ID3) 

4.1 Overview of Decision Trees 

Decision tree learning is a method for approximating discrete-valued target 
functions and is one of the most widely used and practical methods for inductive 
inference. Decision tree learning is generally best suited to problems with the following 
characteristics: Attribute-value pairs represent instances. 

Instances are described by a fixed set of attributes (e.g., temperature) and their 
values (e.g., hot). The easiest situation for decision tree learning occurs when each 
attribute takes on a small number of disjoint possible values (e.g., hot, mild, cold). 
Extensions to the basic algorithm allow handling real-valued attributes as well (e.g., a 
floating point temperature). The target function has discrete output values. A decision 
tree assigns a classification to each example. Simplest case exists when there are only 
two possible classes (Boolean classification). Decision tree methods can also be easily 
extended to learning functions with more than two possible output values. A more 
substantial extension allows learning target functions with real-valued outputs, although 
the application of decision trees in this setting is less common. Disjunctive descriptions 
may be required. Decision trees naturally represent disjunctive expressions. 

4.2 Decision Tree Representation 

A decision tree is an arrangement of tests that prescribes an appropriate test at 
every step in an analysis. In general, decision trees represent a disjunction of 
conjunctions of constraints on the attribute-values of instances. Each path from the tree 
root to a leaf corresponds to a conjunction of attribute tests and the tree itself to a 
disjunction of these conjunctions. More specifically, decision trees classify instances by 
sorting them down the tree from the root node to some leaf node, which provides the 
classification of the instance. Each node in the tree specifies a test of some attribute of 
the instance, and each branch descending from that node corresponds to one of the 
possible values for this attribute. 
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An instance is classified by starting at the root node of the decision tree, testing 
the attribute specified by this node, then moving down the tree branch corresponding to 
the value of the attribute. This process is then repeated at the node on this branch and so 
on until a leaf node is reached. 

4.3 Decision Tree Diagram 

A Each nonleaf node is connected to a test that splits its set of possible 
answers into subsets corresponding to different test results. 

A Each branch carries a particular test result's subset to another node. 

^ Each node is connected to a set of possible answers. 



Set of possible answers Set of possible answers 

Fig 4.1 Tree Structure 

4.4 Quinlan's ID3 Algorithm for Constructing a Decision Tree 
General Form 

Until each leaf node is populated by as homogeneous a sample set as 
possible. Select a leaf node with an inhomogeneous sample set. Replace that leaf 
node by a test node that divides the inhomogeneous sample set into minimally 
inhomogeneous subsets, according to an entropy calculation. 

Specific Form 

Examine the attributes to add at the next level of the tree using an entropy 
calculation. Choose the attribute that minimizes the entropy. 
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4.4.1 Choosing an Attribute using the Entropy Formula 

The central focus of the ID3 algorithm is selecting which attribute 
to test at each node in the tree. 

Procedure: 

See how the attribute distributes the instances. 

Minimize the average entropy. 

Calculate the average entropy of each test attribute and 
choose the one with the lowest degree of entropy 

4.4.2 Entropy Formulae 

Entropy, a measure from information theory, characterizes the (im) purity, 
or homogeneity, of an arbitrary collection of examples. 

Given: 

rib = the number of positive instances in branch b. 

ribc = the total number of instances in branch b of class c. 

n, = the total number of instances in all branches. 

Probability 

P b = Probability an instance on a branch bi$ positive 

number of positive instances on branch %■ 

total number of instances on branch 

If all the instances on the branch are positive, then P b = 1 
(homogeneous positive) If all 

the instances on the branch are negative, then P b = 0 (homogeneous 
negative) 

Entropy 

ip / ■%? ■ %% / ^ 

Entropy = f (— % ? 
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H=I-(P b log 2 (P b )) ■ 

Average Entropy (H) 


•4.1 


Or H(C\A k ) = ^P(a k ,j) 


./=! 


N 


* 


~'ZP<,C l \a kj ylo g2 P(C l \a lJ ) 


■4.2 


Where, 

H (C|A k ) = entropy of classification property of attribute A k . 

P ( ak .) = probability of attribute k being at value j. 

p (q | a| . .) = probability that the class value is Q when attribute k is at its 
j th value. 

]y[ k = Total number of values for attribute A k . j — 1,2, . . ., M k . 

]\f= total number of different classes, i — 1, 2. . . N. 

K= total number of attributes; k -1 , 2. . . K. 

4.5 Use of ID3 in developing rule base for RSP 

Two different sets of rule bases were formed for implementation of the system. 

> Rule Base from Expert Advice: The rule base was formulated on the knowledge 
acquired from engineers and worker at the steel plant right from the supervisors 

downward during visit at the factory premises. 

> Data based: This rule base was formed with the data made available from plant. 
The rules were formulated by assigning discrete values of high, medium and low 
as per the numeric values taken by them. 

4.5.1Rule Base from Expert Advice 

As the available data at the start of project was very few in number and 
insufficient to develop the rule base, the development of expert system was 
started with rule base made on the advice, and experience of the experts. As the 
number of input parameters was high (10) and each of the attributes above took 
values like very high, high, medium, low and very low the number of rules to get 
the green bulk density would be enormous ( 3 8 * 5 2 to be precise as PMT and 
GBD had 5 values and remaining 8 attributes had 3 values). Thus the rule sets 
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were formulated so as to reduce the memory intensive computing that would be 
involved otherwise and the entire process of brick manufacture was broken up 


into sub processes as shown in the table 4.1 below, 


Sub process 1 

Inputs-3 

Output- 1 


CTI 

GT 


GHDT 



RT 


Sub process2 

Inputs-3 

Output- 1 


GT 

MIX TEMP 


MIXER TEMP 



MIXING TIME 


Sub process3 

Inputs-3 

Output- 1 


MIX TEMP 

PMT 


ST 



No of ROLLS 


Sub process4 

Inputs-4 

Output- 1 


PMT 

GBD 


PITCH 



GRAPHITE 



FP 


Sub process5 

Inputs-3 

Output- 1 


GBD 

CCS 


PITCH 



TEMPERING TEMP 


Sub process6 

Inputs-2 

Output- 1 


GBD 

CP 


PITCH 


Sub process7 

Inputs-3 

Output- 1 


CP 

SCR 


Carbon Content 



Chemical Purity 


Sub process8 

Inputs-3 

Output- 1 


GBD 

OXI RESISTANCE 


Carbon Content 



Metal Powder 


Sub process9 

Inputs-3 

Output- 1 


GBD 

HMOR 


Carbon Content 



CCS 



Metal Powder 


Table 4 

1.1 Sub processes for rule generation 
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4.5.2 Sub Process-1 


Let us go on to solve the expert base data set on GT using ID3 algorithm 
and construct the tree structure and generate the rules. 

Entropy calculations 

The data set based on expert advice is shown in table 4.2 of AppendixA in 
which we have input attributes as CTI, GHDT and RT, hence k=3, 

There are 5 classes (veryjhigh, high, medium, low & very_low) for the 
output, hence N=5. 

There are 3 classes each for the inputs (high, medium & low). Mi, M 2 and 
M 3 =3 

Using equation 4.2 we can compute the entropy for each attribute, starting 
with CTI we proceed as follows, 

P (CTI = high) = 9/27 =1/3 
P (CTI = medium) = 9/27 =1/3 
P (CTI = low) = 9/27 =1/3 
P (GT = very high j CTI= high) = 1/9 
P (GT=high | CTI=high) = 2/9 
P (GT=medium | CTI=high) = 4/9 
P (GT=low | CTI=high) = 1/9 
P (GT=very low | CTI=high) = 1/9 
P (GT=very high | CTI=medium) = 1/9 
P (GT=high J CTI=medium) = 2/9 
P (GT=medium | CTI=medium) = 3/9 
P (GT=low| CTI=medium) = 2/9 
P (GT=very low | CTI=medium) = 1/9 
P (GT=very high j CTI=low) = 0/9 
P (GT= high | CTI=low) = 1/9 
P (GT=medium | CTI=low) = 3/9 
P (GT=low | CTI=low) = 3/9 
P (GT=very low | CTI=low) = 2/9 
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Similarly corresponding values for GHDT and RT are computed and we 
have by substituting in the equation 4.2, 

H (GT | CTI) = 0.6168139456208 
H (GT | GHDT) = 0.4873750559284 
H (GT | RT) = 0.5330855183249 

Thus we find that GHDT has the minimum average entropy and hence a 
decision tree with root node as GHDT is constructed as shown in Fig4.2 



H, M, L, VL 


Fig4.2 Incomplete Branching 


Listed below each attribute value in the fig4.2 are the classes for Grain 
temp associated with that particular branching of decision tree, for example under 
GHDT equals High we have all classes for GT where as under GHDT equals Low 
we have all classes except very high. Since all three branches result in an 
inconclusive classification, we have to branch in yet another attribute. 

Now we must partition table 4.2 according to attribute GHDT. This results 
in three subtables (one each for GHDT equals high, medium and low). The 
subtable for GHDT equals high is shown in Table (4.2A). Once again we compute 
the entropy for each attribute in the three subtables. The other two sub tables are 
shown in Table (4.2 B & C). 
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CT1 

RT 

GRAIN TEMP 

HIGH 

HIGH 

VERY HIGH 

HIGH 

MEDIUM 

HIGH 

HIGH 

LOW 

MEDIUM 

MEDIUM 

HIGH 

VERY HIGH 

MEDIUM 

MEDIUM 

HIGH 

MEDIUM 

LOW 

MEDIUM 

LOW 

HIGH 

HIGH 

LOW 

MEDIUM 

MEDIUM 

LOW 

LOW 

LOW 


Table (4.2A) Classification when GHDT equals High 


CTI 


HIGH 

HIGH 

HIGH 

HIGH 

MEDIUM 

MEDIUM 

HIGH 

LOW 

MEDIUM 

MEDIUM 

HIGH 

HIGH 

MEDIUM 

MEDIUM 

MEDIUM 

MEDIUM 

LOW 

LOW 

LOW 

HIGH 

MEDIUM 

LOW 

MEDIUM 

MEDIUM 

LOW 

LOW 

LOW 


Table (4.2B) Classification when GHDT equals Medium 



lillliiplliill! 

HIGH 

HIGH 

MEDIUM 

HIGH 

MEDIUM 

LOW 

HIGH 

LOW 

VERY LOW 

MEDIUM 

HIGH 

MEDIUM 

MEDIUM 

MEDIUM 

LOW 

MEDIUM 

LOW 

VERY LOW 

LOW 

HIGH 

LOW 

LOW 

MEDIUM 

VERY LOW 

LOW 

LOW 

VERY LOW 


Table (4.2B) Classification when GHDT equals Low 
Entropy Calculations GHDT (High) 

Using equation 4.2 we get, 

H(GT| CTI)=3/9{ 3 * [ 1 /3 * lo g( 1 /3 )]+ [ 1 /3 * log( 1 /3 )+2/3 log(2/3 )] 

+ 3* [l/3*log(l/3)]}= 0.477121254719662437295 
H(GT| RT)=3/9{ 3*[l/3*log(l/3)]+ 3*[l/3*log(l/3)] 

+[l/3*log(l/3)+2/31og(2/3)]}= 0.2764345909436749738 
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Here we see the RT has the lowest entropy. Hence it will form the next node and 
since only one attribute is left that of CTI which will form the terminal node and the tree 
structure will be complete as in fig 4.3. and fig 4.4. 



Fig4.3 Partially complete Branching 
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Entropy Calculations GHDT (Medium) 

Using equation 4.2 vve get, 

H(GT| Cri)=3/9{ 3 * [ 1 /3 *log( 1 /3)]+[ 1 /3 * log( 1 /3)+2/3 log(2/3)] 

+ 3 ,fc [l/3*log(l/3)]}= 0.34333014553567079497 
H(GT| RT)=3/9{ 3*[l/3*log(l/3)]+3*[l/3*log(l/3)] 

+[ 1/3 * log ( 1 /3)+2/3 log(2/3 )] } = 0.18428972729578331587 

Entropy Calculations GHDT (Medium) 

Using equation 4 .2 we get, 

H(GT| CTI)=3/9{ 3*[l/3*log(l/3)]+[l/3*log(l/3 )+2/3 log(2/3 )] 

4- 3*[l/3*log(l/3)]}= 0.41022570012766661613 
H(GT| RT)=3/9{ 3*[l/3*log(l/3)]+ 3*[l/3*log(l/3)] 

+[l/3*log(l/3)+2/31og(2/3)]}= 0.18428972729578331587 
Thus we see that for GHDT Medium and Low also we get RT having lower 
entropy than CTI. Software has been used to calculate the average entropy using 
Quinlan's ID3 Algorithm and to generate the decision trees for the data for the processes 
attached as Table 4.3 to 4.10 in AppendixA. The trees generated for data sets from Sub 
processes 1 to 9 has been attached as fig 4.5 to 4.13 in Appendix B. 

4.5.3 Generation of Rule Base from Available Data 

The Data available at the start from plant was only 45 complete sets up to the 
Green BD stage. Subsequently further data was provided up to 86 sets. The Data Based 
Rules have been formed only up to the green BD stage as at present no data is available 
for Sub processes beyond GBD stage. The rules were made with the following inputs and 
Green BD as the output: 

• Combustion Chamber Temperature 

• Grain Heater Drum Temperature 

• Retention T ime in Grain Heater 

• Mixer temperature 

• Mixing time 

• Mix temperature 

• Storage time 

• Number of rolls 

• Press Mix T emperature 

• Forming pressure 
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The other parameters like the granulometry, pitch quantity and graphite quantity 
in the data was fixed and remained constant. Thus the same were not used as input. The 
86 data sets shown in table 4. 1 1 A of Appendix A was converted to a table of very high, 
high, medium, low and very low as the case may be by taking calculating the mean and 
standard deviation(SD) of the attribute values, after which values below 2 SD were 
assigned to low, within 2 SD to Medium and all above 2 SD to High, and the same is 
shown in table 4.1 IB of Appendix A. The IDS algorithm was applied to this data set in 
table 4. 1 IB. The tree structure obtained for the data is shown in fig 4. 14of Appendix B. 

The results obtained from the data was in variation with the expert opinion 
because the data collected had very high noise content and secondly the quantity of data 
sets was far too small to capture major portion of the feature space. It was later admitted 
by the plant staff that the data could be noisy due lack of proper instrumentation and 
ignorance on part of the ground workers on the importance of collecting accurate data. To 
observe a pattern with in the data at least 3000 to 4000 sets of data will be required. The 
system can be made robust and reliable only once it is able to capture 90 to 95 % of the 
feature space. For the 86 data sets the IDS algorithm generated 66 sets of rules that can be 
read out from the tree structure attached as per fig 4.14 of Appendix B. 
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CHAPTER 5 

CLASSIFICATION AND REGRESSION TREES 


5.1 Overview of cart 

CART is a robust data-mining and data-analysis tool that automatically searches for 
important patterns and relationships and quickly uncovers hidden structure even in highly 
complex data. A decision tree is a flow chart or diagram representing a classification system or 
predictive model. The tree is structured as a sequence of simple questions, and the answers to 
these questions trace a path down the tree. The end point reached determines the classification or 
prediction made by the model. 

The CART methodology is technically known as binary recursive partitioning. The 
process is binary because parent nodes are always split into exactly two child nodes and 
recursive because the process can be repeated by treating each child node as a parent. The key 
elements of a CART analysis are a set of rules for: 

• Splitting each node in a tree, 

• Deciding when a tree is complete, and 

• Assigning each terminal node to a class outcome (or predicted value for 
regression). 

5.2 Splitting Rules 

To split a node into two child nodes, CART always asks questions that have a “yes” or 
“no” answer. For example, in the tree developed for RSP using the 117 data sample shown in fig 
5.1 below. 



Fig 5.1 Tree structure developed for RSP using 117 data sets. 
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Node 1 was split on ST, a case goes left if ST <= 58.500, i.e. the root, or parent node is 
split into two branches with “yes” cases going to the left child node and “no” cases to the right, 
and Node 2 was split on GHDT, a case goes left if GHDT <= 537.500, similarly Node 3 was 
split on CTI, a case goes left if CTI <= 342.000. 

5.3 Choosing a Split 

Next step is to rank order each splitting rule on the basis of a quality-of split criterion. 
The default criterion used in CART is the GINI rule; essentially a measure of how well the 
splitting rule separates the classes contained in the parent node. Five alternative criteria are also 
available for classification trees and two criteria for regression trees. In addition, to deal more 
effectively with select data patterns, CART also offers splits on linear combination of continuous 
predictor variables. For example shown in Fig 5.2 below is partial portion of the tree developed 
up to level 3 in the case of RSP data the best quality slip obtained at Node 1 slips all 117 data, 
Node 2 slips 82 data sets , Node 3 splits 62 sets is done using GINI rule. 



5.4 Class Assignment 

Once a best split is found, CART repeats the search process for each child node, 
continuing recursively until further splitting is impossible or stopped. Splitting is impossible if 
only one case remains in a particular node or if all the cases in that node are exact copies of each 
other (on predictor variables). CART also allows splitting to be stopped for several other reasons, 
including that a node has too few cases. (The default for this lower limit is 10 cases, but may be 
set higher or lower to suit a particular analysis). Once a terminal node is found we must decide 
how to classify all cases falling within it. One simple criterion is the plurality rule: the group 
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Classification tree topology for: GBD 


with the greatest representation determines the class assignment. For example shown in Fig 5.3 
below is the case for class 1 or very low GBD in red square boxes, the degree of redness gives 
the class assignment for it. The best to worst case changing gradually from red to blue, in the fig 
Below only three nodes have cases of very low GBD. 



Fig 5.3 Class Assignment based on Number of cases in a Node 
5.5 Testing 

Once the maximal tree is grown and a set of sub-trees is derived from it, CART 
determines the best tree by testing for error rates or costs. With sufficient data, the simplest 
method is to divide the sample into learning and test sub-samples. The learning sample is used to 
grow an overly large tree. The test sample is then used to estimate the rate at which cases are 
misclassified (possibly adjusted by misclassification costs). The misclassification error rate is 
calculated for the largest tree and also for every sub-tree. The best sub-tree is the one with the 
lowest or near-lowest cost, which may be a relatively small tree. Some studies, as the one in our 
case of RSP data, where there are insufficient data samples to allow a good-sized separate test 
sample, CART employs the computer-intensive technique of cross validation. 
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5.6 Cross Validation 

In such cases, CART grows a maximal tree on the entire learning sample. This is the tree 
that will be pruned back. CART then proceeds by dividing the learning sample into 1 0 roughly - 
equal parts, each containing a similar distribution for the dependent variable, in the case of 117 
data sets each sample would have 12 samples with the last having 8 samples. CART takes the 
first 9 parts of the data, constructs the largest possible tree, and uses the remaining 1/10 of the 
data to obtain initial estimates of the error rate of selected sub-trees. The same process is then 
repeated (growing the largest possible tree) on another 9/10 of the data while using a different 
1/10 part as the test sample. The process continues until each part of the data has been held in 
reserve one time as a test sample. The results of the 10 mini-test samples are then combined to 
form error rates for trees of each possible size; these error rates are applied to the tree based on 
the entire learning sample. 

The upshot of this complex process is a set of fairly reliable estimates of the independent 
predictive accuracy of the tree. This means that we can know how well any tree will perform on 
completely fresh data — even if we do not have an independent test sample. Because the 
conventional methods of assessing tree accuracy can be wildly optimistic, cross validation is the 
method CART normally uses to obtain objective measures for smaller data sets. 

5.7 Surrogate Splitter 

CART handles missing values in the database by substituting “surrogate splitters,” which 
are back-up rules that closely mimic the action of primary splitting rules. Suppose that, in a given 
model, CART splits data according to household income. If a value for income is not available, 
CART might substitute education level as a good surrogate. 

The surrogate splitter contains information that is typically similar to what would be 
found in the primary splitter. Other products’ approaches treat all records with missing values as 
if the records all had the same unknown value; with that approach all such “missing” are 
assigned to the same bin. In CART, each record is processed using data specific to that record; 
this allows records with different data patterns to be handled differently, which results in a better 
characterization of the data. 
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5.8 Special Features of CART 

Easy Data Access 

With a direct link to DBMS, CART provides easy access to over 80 different file 
formats and supports ODBC for both reading and writing. For example, you can 
import and export statistical analysis packages (e.g., SAS, SPSS), databases (e.g., 
Oracle, Informix), and spreadsheets (e.g., Excel, Lotus). 

Automatic Rule Generation 

With a click of a button one can, view if-then rules for any individual node 
or for the entire tree. To facilitate external classifying and scoring of new data, the 
splitting rules are formatted as C-compatible programming statements that can be 
exported as either a text file or cut and pasted into another application. 

5.9 Generation of Cart Tree for RSP data. 

The opening screen of the cart programme displays as under, which is very similar to the 
windows icons, the function of each of the buttons is explained below 



Fig 5.4 the opening screen of Cart. 

The most commonly-used commands have corresponding toolbar icons. The following 
icons are most frequently used: 

Opening an existing file for analysis 

Submitting a text file of CART commands for batch processing 
Turning command-line entry mode on or off 
Opening a log of processed CART commands 
Setting reporting, random number and directory options 
Specifying model variables and parameters that control analysis 
Starting CART analysis of current model 
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When the open button is pressed one is prompted to load the data file, in which the RSP 


data set in Excel format is loaded to give the following window in Fig 5.5 below, 



Advancet 


1 Penalty J 

Select Cases Y I est T ree 


Model 


Categorical 
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Select Target and Predictor Variables 
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MIXERJTE 
M IXI N G_T 
MIXJEMP 
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C 4 Regression 


Weight Variable 


Save T ree Information. 


Cancel 


Fig 5.5 Loading RSP data into Cart 



Fig 5.6 Tree generated for RSP data 


48 








The target variable selected is GBD and the rest are the predictors. The number of levels 
is set to 5 for the target variable corresponding to 1 to 5 being very low, low, medium, high and 
very high, hence minimum value is set to 1 . Then the grow tree button is pressed to give the 
following screen as in Fig 5.6 above. 

5.10 Results obtained and use of Cart analysis in RSP data 

Since the data collected from the plant was noisy in nature and not in adequate quantity 
hence the results obtained did not completely match with the expert opinion, which was 
adequately acknowledged by the RSP staff also. The salient features are discussed below. 

The rule generated by cart is attached as per Appendix D. The rules generated are only 21 
in number for the given data as compared to 66 rules generated by ID3, and 243 rules generated 
by matlab, for the same amount of data. However the ranges for high, medium, low are in 
general agreement with that given to us by the experts with minor changes. 

The variable importance as deduced by Cart was not in agreement with the expert opinion 
due to reasons just brought earlier. The comparison of variable importance and the other out puts 
obtained by Cart Analysis is attached as Fig 5.7 & 5.8 of Appendix E. 
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CHAPTER 6 


BAYESIAN NETWORKS AND PROBABILTY PREDICTIONS 

6.1 Overview of Bayesian networks 

The essence of the Bayesian approach is to provide a mathematical rule 
explaining how you should change your existing beliefs in the light of new 
evidence. In other words, it allows us to combine new data with our existing 
knowledge or expertise. Mathematically, 

Bayes’ rule states 


Posterior _ Pr obabilty 


Conditional _ Pr obabilty * Aprior _ Pr obablity 
Pr obability 


rs ■ u , run , n P(e \ R = r)P(R = r) 

Or, in symbols, P(R = r | e) = ~P() 

Where P(R=r|e) denotes the probability that random variable R has value r 
given evidence e. 

The denominator is just a normalizing constant that ensures the posterior 
adds up to 1; it can be computed by summing up the numerator over all 
possible values of R, i.e., 

P (e) = P(R=0, e) + P(R=1, e) + ... = sum P(e | R=r) P(R=r) 

For complicated probabilistic models, computing the normalizing constant P (e) is 
computationally intractable (because there are an exponential or even infinite 
cumber of values of “R” to sum over); graphical models help simplify these 
.putations so that they become feasible. 


Graphical Models 


Graphical models are a marriage between probability theory and graph theory, 
y provide a natural tool for dealing with two problems that occur throughout 
lied mathematics and engineering — uncertainty and complexity — and in 
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particular they are playing an increasingly important role in the design and analysis 
of machine learning algorithms. Fundamental to the idea of a graphical model is the 
notion of modularity — a complex system is built by combining simpler parts. 
Probability theory provides the glue whereby the parts are combined, ensuring that 
the system as a whole is consistent, and providing ways to interface models to data. 
The graph theoretic side of graphical models provides both an intuitively appealing 
interface by which humans can model highly interacting sets of variables as well as 
a data structure that lends itself naturally to the design of efficient general-purpose 
algorithms. 

6.2.1 Types of Graphical Models 

Graphical models are graphs in which nodes represent random variables, 
and the lack of arcs represents conditional independence assumptions. 
Undirected graphical models, also called Markov Random Fields 
(MRFs), have a simple definition of independence: two (sets of) nodes A 
and B are conditionally independent given a third set, C, if all paths between 
the nodes in A and B are separated by a node in C. By contrast, directed 
graphical models, also called Bayesian Networks or Belief Networks 
(BNs), have a more complicated notion of independence, which takes into 
account the directionality of the arcs. (Note that, despite the name, Bayesian 
networks do not necessarily imply a commitment to Bayesian methods; 
rather, they are so called because they use Bayes' rule for inference). 
Undirected graphical models are more popular with the physics and vision 
communities, and directed models are more popular with the AI and 
statistics communities. 

6.3 Directed Graphical Models 
6.3.1 Bayesian Networks 

Although directed models have a more complicated notion of 
independence than undirected models, they do have several advantages. The 
most important is that one can regard an arc from A to B as indicating that A 
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causes B. This can be used as a guide to construct the graph structure. In 
addition, directed models can encode deterministic relationships, and are 
easier to learn (fit to data). 



Fig 6.1 Graph Structure and Conditional Probability Table 

(CPT) 


In addition to the graph structure, it is necessary to specify the 
parameters of the model. For a directed model, we must specify the 
Conditional Probability Distribution (CPD) at each node. If the variables are 
discrete, this can be represented as a Conditional Probability table (CPT), 
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which lists the probability that the child node takes on each of its different 
values for each combination of values of its parents. Consider the following 
example in the Fig 6.1 above, in which all nodes are binary, i.e., have two 
possible values, which we will denote by T (true) and F (false). We see that 
the event "grass is wet" (W=true) has two possible causes: either the water 
sprinker is on (S=true) or it is raining (R=true). For example, we see that Pr 
(W=true | S=true, R=false) = 0.9 (second row), and hence, Pr (W=false | 
S=true, R=lalse) = 1 - 0.9 = 0.1, since each row must sum to one. Since the 
C node has no parents, its CPT specifies the prior probability that it is 
cloudy (in this case, 0.5). 

6.3. 1.1 Conditional Independence 

The simplest conditional independence relationship encoded in a 
Bayesian network can be stated as follows: a node is independent of its 
ancestors given its parents, where the ancestor/parent relationship is with 
respect to some fixed topological ordering of the nodes. By the chain rule of 
probability, the joint probability of all the nodes in the graph above is P(C, 
S, R, W) = P (C) * P (S|C) * P (R|C, S) * P (W|C, S, R) 

by using conditional independence relationships, we can rewrite this 
as P(C, S, R, W) = P(C) * P (S|C) * P (R|C) * P (W|S, R). Where we 
were allowed to simplify the third term because R is independent of S given 
its parent C-, and the last term because W is independent of C given its 
parents S and R. We can see that the conditional independence relationships 
allow us to represent the joint more compactly. Here the savings are 
minimal, but in general, if we had n binary nodes, the full joint would 
require 2" space to represent, but the factored form would require n* 2 k ) 
space to represent, where k is the maximum fan-in of a node. And hence 
fewer parameters make learning easier. 

6.3. 1.2 Inference 

The most common task we wish to solve using Bayesian networks is 
probabilistic inference. For example, consider the water sprinkler network, 
and suppose we observe the fact that the grass is wet. There are two possible 
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likely? We can use Bayes' rule to compute the posterior probability of each 
explanation (where 0==false and l==true). 


P(S = 1 | W = 1) = = 

P(W = 1) 

I ,,-P(C = c,S - H* = r,W = 1) 0.2781 

P(W = 1) ~ 0.6471 


0.430 


P(R = 1 | W = 1) = 


pjj. = ijr = i) 
P(W = 1) 
s,R = 1 ,W = 1) 


P{W = 1) 


Where, 


0.4851 

0.6471 


0.708 


P(W = 1) = £ F ( C = c,S = s, R - r,W = 1) = 0.6471 

c ,r , s 


is a normalizing constant, equal to the probability (likelihood) of the 
data. So we see that it is more likely that the grass is wet because it is 
raining: the likelihood ratio is 0.7079/0.4298 = 1.647. 

In the above example, notice that the two causes “competes" 
to "explain" the observed data. Hence S and R become conditionally 
dependent given that their common child, W, is observed, even 
though they are marginally independent. For example, suppose the 
grass is wet, but that we also know that it is raining. Then the 
posterior probability that the sprinkler is on goes down: P (S=1|W=T, 
R=l) = 0.1945, this is called "explaining away". 

6.3.1.3 Top-Down and Bottom-Up Reasoning 

In the water sprinkler example, we had evidence of an effect 
(wet grass), and inferred the most likely cause. This is called 
diagnostic, or "bottom up", reasoning, since it goes from effects to 
causes; it is a common task in expert systems. Bayes nets can also be 
used for causal, or "top down", reasoning. For example, we can 
compute the probability that the grass will be wet given that it is 
cloudy. Hence Bayes nets are often called "generative" models, 
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because they specify how causes generate effects. The particular 
aspect also termed as backward chaining can be exploited for setting 
the operating parameters as the desired out put is well known, thus in 
order to “effect” the desired out put the Bayesian tree can be used to 
find out the causes, this aspect can also be used for trouble shooting . 
In case the output obtained is not as per the desired specification one 
can find out the “cause”. 

6.4 Generation of Tree from RSP data. 

In order to generate the tree the javabayes software has been used. When the 
software is loaded and run it generates two windows the javabayes editor and the 
javabayes console shown in Fig 6.2 below. In the console window is displayed a 
series of instructions as one proceeds to generate the network on the editor window. 
The opening instructions on the javabayes console are as follows, to start editing 
networks, press the Create button and click on the JavaBayes editor, or load a 
network using the file -> Open menu. Thus it is highly user friendly and fully menu 
driven. The editor 

window has the following buttons for editing of the networks, create, move, delete, 
edit variable, edit function and edit networks. The create button is used to create the 
nodes, shown as green circles in the Fig 6.3 below, for this one has to just click the 
create button and then click any where inside the editor window, where one wants 
to create a node. Once all the nodes have been created as desired, the input nodes 
are connected to the output node by drawing arrows (arcs) from input to output 
using the mouse (just click and drag). Once the entire tree model is in place, if one 
wants to move or delete a node or arc, the same can be achieved using the move or 
delete buttons, for this one has to first click the move/delete button and then 
drag/ click the desired node/arc which is to be moved/deleted. 
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JavaBayes Editor 
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the Network->Open menu. 
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Fig6.2 JavaBayes Editor and Java Console 
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Fig 6.3 Tree Structure for GBP & Labeling Parameters for GHDT 
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Once the graph structure has been finalized, it is necessary to specify the parameters 
for the model. We must specify the Conditional Probability Distribution (CPD) at 
each node. This can be represented as a Conditional Probability table (CPT), which 
lists the probability that the child node takes on each of its different values for each 
combination of values of its parents. The operation is achieved using the edit 
variable and edit function buttons. Here again the edit variable button is clicked and 
then the variable to be edited is clicked which throws up a window as shown in the 
right side of Fig 6.3. From the window on the right side it can be seen that it allows 
you to change the name of the variable and the range of values that the variable can 
take as high, medium, low etc. 

Once this is done for all the nodes of the tree one has to start entering the 
values for the Conditional Probability table (CPT). This done using the edit function 
button, when any node is clicked after clicking the edit function button it throws up 
a window as shown in Fig 6.4 below, which is for the node Grain Temperature. 


[|||§Edit Function 
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Fig 6.4 Edit Function Window for Entering the Conditional Probability Table. 


The values entered here are those obtained from the rules generated by the 
experts, which is attached as per table 4.2 to 4.10 of Appendix A. Thus the tree has 
been developed only for the expert based opinion as the volume of data generation 


57 





been developed only for the expert based opinion as the volume of data generation 
for calculation of the probabilities was inadequate to develop a robust and 
meaningful system. To start tasting we have to use the query and observe buttons. 
The moment any of the nodes is clicked after clicking the query button, the 
probability distribution table of the nodes values gets listed out in the console 
window. Now if we observe that 


Probability (“Combusion_chamber_temp”) =1.0// P (Medium | evidence) 
Probability (“Grain_heater_drum_temp”) = 1 .0// P (High | evidence) 
Probability (“Retention time”) =1 .Oil P (Low | evidence); 

When we query for the node grain temperature we get, 

Probability (“Grain_temp”) = 1.0// P (Medium |evidence) shown in Fig 6.5 
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To observe a node, click on it. 

To query on a particular node, click on it. 

Posterior distribution: 

probability ( "Combusion_chamber_temp" ) { //I variable(s) and 3 v 
table 

0.0 // p(High | evidence ) 

1.0 // p(Medium | evidence ) 

0.0; //p(Low | evidence); 

} 

Posterior distribution: 

probability ( ,, Grain_heater_drum_temp ,, ){//l variable(s) and 3 vah 
table 

1.0 // p(High | evidence ) 

0.0 // p(Medium | evidence ) 

0.0; // p(Low | evidence ); 

} 

Posterior distribution: 

probability ( "Retentionjime" ) { //I variable(s) and 3 values 
table 

0.0 // p(High | evidence ) 

0.0 // p(Medium | evidence ) 

1.0; // p(Low 1 evidence ); 

} 

Posterior distribution: 
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Fig 6.5 Probability Prediction for Grain Temp (Medium=l) Given RT, GHDT 


&CTI 


6.5 Application of the Results Obtained for RSP 

We proceed in this manner till we reach GBD, using “observe” and “query” 
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buttons to obtain results as compiled in Fig 6.6 to Fig 6.8 of Appendix F for a 
specific arbitrary example. For the processes following GBD an example has been 
dealt with and compiled as Fig 6.9 to 6.12 of the same Appendix F. Th.us it can be 
seen that the effect of change in parameter in any of the parent node can be 
immediately queried at any level of the child node. Further, tree models have also 
been developed for backward chaining processes to diagnose the faults in case the 
out put shifts from the desired specification limits. 

These Tree models developed for RSP can be used for the following 
purposes in the plant; 

> It can be used for finding the effect of change in any single parameter on the 
final out put for example if want to know what the GBD is likely to be if GHDT 
alone is High, we just have to observe GHDT as High and query for GBD to get 
the result; Posterior distribution: probability ("Green_bd") {1 variable and 5 
values table} 

0.03566529492455419 P (Veryjiigh | evidence) 

0.13397347965249204 P (High | evidence) 

0.293278463648834 P (Medium j evidence) 

0.28321902149062644 P (Low | evidence) 

0.2538637402834934; P (Very_low | evidence); Hence the probability of the 
GBD being Medium or low is highest compared to the other three. 

> Thus it can be used to refine the shape of the membership functions of the 
expert based system to match closely to the realistic data when ever the same is 
made available by the authorities for validation purposes. 

> Validation of the rules from time to time as they are updated with additional 
data is most easily done using these tree structures as it is simple and easy to use 
even by the plant operator. 

> During the validation of data the, these tree structures can he of great help to 
give an overall view of how the system is responding to the inputs. 

During the manufacturing process it can be readily used to predict the out put 
and hence take corrective action immediately. 
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> The trees developed for backward chaining can be used for trouble shooting/ 
fault diagnosis incase of wide variation in the output parameters. 
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CHAPTER 7 

DEVELPOMENT OF FIS FROM EXPERT DATA 


can call from the command 

of 


7.1 Introduction to Fuzzy Toolbox of Matlab 

The Fuzzy Logic Toolbox is a collection of functions built on the MATLAB nume 
computing environment. It provides tools for you to create and edit fuzzy inference system 
within the framework of MATLAB. This toolbox relies heavily on graphical user interface (G 
tools to help you accomplish your work, although you can work entirely from the command h 
if you prefer. 

The toolbox provides two categories of tools: 

•Command line functions 
•Graphical, interactive tools 

The first category of tools is made up of functions that you 

line or from your own applications. Many of these functions are MATLAB M files, se 

MATLAB statements that implement specialized fuzzy logic algorithms. You can vi 

• Tne You can change the 

MATLAB code for these functions using the statement type functionjram • 

way any toolbox function works by copying and renaming the M-file, then modifying yo 

You can also extend the toolbox by adding your own M-files. 

i that let you access many of 

Secondly, the toolbox provides a number of interactive tools uw- 

t for fuzzy 

the functions through a GUI. Together, the GUI- based tools provide an environm n 
inference system design, analysis, and implementation. 

7.2 Building Systems with the Fuzzy Logic Toolbox 

Although it’s possible to use the Fuzzy Logic Toolbox by working strictly 

, . iiv There are five primary 

command line, in general it’s much easier to build a system graphically- Logic 

GUI tools for building, editing, and observing fuzzy inference systems 

' . rshiD Function Editor, the 

Toolbox: the Fuzzy Inference System or FIS Editor, the Membersn p _ 

Rule Editor, the Rule Viewer, and the Surface Viewer. These GUIs are dyn ^ 

that changes you make to the FIS using one of them, can affect what you see y 

open GUIs. You can have any or all of them open for any given system. 

„ . , , AC the graphical ANFIS Editor 

In addition to these five primary GUIs, the toolbox includes me b y 

. ntive neural fuzzy inference 

GUI, which is used for building and analyzing Sugeno-type aaapt 
systems. 
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The FIS Editor handles the high level issues for the system: How many input and output 
variables? What are their names? The Fuzzy Logic Toolbox doesn’t limit the number of inputs. 
However, the number of inputs may be limited by the available memory of your machine. If the 
number of inputs is too large, or the number of membership functions is too big, then it may also 
be difficult to analyze the FIS using the other GUI tools. 

The Membership Function Editor is used to define the shapes of all the membership 
functions associated with each variable. 

The Rule Editor is for editing the list of rules that defines the behavior of the system. 

The Rule Viewer and the Surface Viewer are used for looking at, as opposed to editing, 
the FIS. They are strictly read-only tools. The Rule Viewer is a MATLAB-based display of the 
fuzzy inference diagram shown at the end of the last section of chapter 3. Used as a diagnostic, it 
can show (for example) which rules are active, or how individual membership function shapes 
are influencing the results. The Surface Viewer is used to display the dependency of one of the 
outputs on any one or two of the inputs — that is, it generates and plots an output surface map for 
the system. 

7.3 FIS Developed For RSP 

Two Fuzzy Inference Systems were developed for RSP, one based on the expert opinions 
where all rules were framed by taking inputs from the engineers and works at RSP and RDCIS, 
and the other was based on the data collected at the plant. In this chapter the former system will 
be discussed and the later will be discussed in chapter 8. The Expert based system was developed 
using fuzzy editor, while the data based system was developed using ANFIS editor. 

The fuzzy inference engine developed for RSP using the tools mentioned in section 7.2 is 
shown in Fig 7.1 below. The five primary GUIs can all interact and exchange information. Any 
one of them can read and write both to the workspace and to the disk (the read-only viewers can 
still exchange plots with the workspace and/or the disk). For any fuzzy inference system, any or 
all of these five GUIs may be open. If more than one of these editors is open for a single system, 
the various GUI windows are aware of the existence of the others, and will, if necessary, update 
related windows. Thus if the names of the membership functions are changed using the 
Membership Function Editor, those changes are reflected in the rules shown in the Rule Editor. 
The editors for any number of different FIS systems may be open simultaneously. The FIS 
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Editor, the Membership Function Editor, and the Rule Editor can all read and modify the FIS 
data, but the Rule Viewer and the Surface Viewer do not modify the FIS data in any way. 

FIS EDITOR 





RULE EDITOR 


MEMBERSHIP FN EDITOR 


Fuzzy 

Inference 

System 


READ ONLY TOOLS 


RULE VIEWER 


SURFACE VIEWER 




fSci 3 |sf" 


j 


Fig 7. 1 Fuzzy Inference System for RSP Showing the Five Tools In the GUI 





For the purpose of developing the FIS model for RSP, the model of 9 sub 
processes mentioned in chapter4 table 4.1 was used. 

7.3.1 Sub process 1 

I shall describe this process with 3 inputs and one out put in detail. In this sub 
process the ranges for the input and output as given by the experts is given in table 7.1 


Parameter 

Range 

CTI 

150-650 degree centigrade 

GHDT 

200-800 degree centigrade 

RT 

2-15 minutes 

GRAIN TEMPERATURE 

90-300 degree centigrade 


Table 7.1 Ranges for Parameters of Sub Processes - 1 


The rules for the system are given in table 4.2 of Appendix A, the ranges for high, 
medium and low have been arrived at after a lot of deliberations between the various 
experts so as to obtain near true data results. Now we need to know the variation of the 
output, the shape of the curves, with respect to the three inputs, by plotting the curve for 
CTI Vs GT, GHDT Vs GT AND RT Vs GT we get the following plots given in Figures 
7.2, 7.3 and 7.4. The remaining set of data are attached as Table 4.3 to 4.10 to Appendix 
A 
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Fig 7.4 Plots for RT Vs GT 

Now that we know the rules, and we know how the output varies with the inputs, 
let’s start working with the GUI tools to construct a fuzzy inference system for the 
decision process of finding what the Grain Temperature is given what the three inputs 
are. 

7.3.2 The Fis Editor 

The FIS Editor displays general information about a fuzzy inference system. 
There’s a simple diagram at the top that shows the names of each input variable on the 
left, and those of each output variable on the right. The sample membership functions 
shown in the boxes are just icons and do not depict the actual shapes of the membership 
functions. 

Below the diagram is the name of the system and the type of inference used. The 
default, Mamdani-type inference, is what we’ve been describing so far. Another slightly 
different type of inference, called Sugeno-type inference, is also available. This method 
has already been explained in “Sugeno-Type Fuzzy Inference” in chapter j and which 
has been used in the ANFIS system to develop the data based FIS system explained. 


-jg|x| 



66 



In the fig 7.5 just below the name of the fuzzy inference system, on the left side 
of the figure, are the pop-up menus that allow you to modify the various pieces of the 
inference process. On the right side at the bottom of the figure is the area that displays the 
name of an input or output variable, its associated membership function type, and its 
range. The latter two fields are specified only after the membership functions have been. 
Below that region are the Help and Close buttons that call up online help and close the 
window, respectively. At the bottom is a status line that relays information about the 
system. 

To start this system from scratch, type fuzzy at the MATLAB prompt. The 
generic untitled FIS Editor opens, with one input, labeled inputl, and one output, labeled 
outputl. For this sub system, we will construct a three-input, one output system, so go to 
the Edit menu and select Add input twice. A second and third yellow box labeled input2 
and input3 will appear. The three inputs we will have in our example are service and 
food. Our one output is tip. We’d like to change the variable names to reflect that, 
though: 

> Click once on the left-hand (yellow) box marked inputl (the box will be 
highlighted in red). 

> In the white edit field on the right, change inputl to CTI and press 

Return. 
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> dick once on the left-hand (yellow) box marked input2 (the box will be 
highlighted in red). 

> In the white edit field on the right, change input2 to GHDT and similarly 
input3 to RT and press Return. 

> Click once on the right-hand (blue) box marked outputl. 

> In the white edit field on the right, change outputl to GT. 

7.3.3 Membership Function Editor 

Next we have to define the membership functions associated with each of the 
variables. To do this, open the Membership Function Editor. You can open the 
Membership Function Editor in one of three ways: 

• Pull down the View menu item and select Edit Membership Functions..., 

• Double-click on the icon for the output variable, GT. 

• Type mfedit at the command line. 


Those menu items allow you 
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The Membership Function Editor in fig 7.6 shares some features with the FIS 
Editor. In fact, all of the five basic GUI tools have similar menu options, status lines, and 
Help and Close buttons. The Membership Function Editor is the tool that lets you display 
and edits all of the membership functions associated with all of the input and output 
variables for the entire fuzzy inference system. 

On the upper left side of the graph area in the Membership Function Editor is a 
“Variable Palette” that lets you set the membership functions for a given variable. To set 
up your membership functions associated with an input or an output variable for the FIS, 
select an FIS variable in this region by clicking on it. 

Next select the Edit pull-down menu, and choose Add MFs. A new window will 
appear which allows you to select both the membership function type and the number of 
membership functions associated with the selected variable. The built in membership 
functions available in the tool box are triangular, trapezoidal, gbell, Gaussian, sigma, 
dsigma, psigma, pi, S and Z membership functions. In the lower right corner of the 
window are the controls that let you change the name, type, and parameters (shape), of 
the membership function, once it has been selected. The membership functions from the 
current variable are displayed in the main graph. These membership functions can be 
manipulated in two ways. You can first use the mouse to select a particular membership 
function associated with a given variable quality, (such as high, for the variable, CTI), 
and then drag the membership function from side to side. This will affect the 
mathematical description of the quality associated with that membership function for a 
given variable. The selected membership function can also be tagged for dilation or 
contraction by clicking on the small square drag points on the membership function, and 
then dragging the function with the mouse toward the outside, for dilation, or toward the 
inside, for contraction. This will change the parameters associated with that membership 
function. 

Below the Variable Palette is some information about the type and name of the 
current variable. There is a text field in this region that lets you change the limits of the 
current variable’s range (universe of discourse) and another that lets you set the limits of 
the current plot (which has no real effect on the system). The process of specifying the 
input membership functions for this two input tipper problem is as follows: 
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> Select the in P ut variable, CTI, by double-clicking on it. Set both the Range and 
the Display Range to the vector [150 650]. 

> Select Add MFs. from the Edit menu. The window below pops open. 

> Use the pull-down tab to choose trimf for MF Type and 3 for Number of MFs. 
This adds three triangle MFs to the input variable CTI. 

'r Click once on the leftmost triangle. Change the name of the curve to low. To 
adjust the shape of the membership function, either use the mouse, as described 
above, or type in a desired parameter change, and then click on the membership 
function. The parameter listing for this curve is [300 350] 

> Name middle triangle, medium, and the rightmost triangle, high. Reset the 
associated parameters to [300 350 400] and [350 400]. 

> Select the input variable, GHDT, by clicking on it. Set both the Range and the 
Display Range to the vector [200 800]. 

> Select Add MFs. from the Edit menu and add three trimf curves to the input 
variable GHDT. 

> Similarly name the triangles as low, medium and high with parameters [300 400], 
[300 400 500] and [400 500] 

> Select the input variable, RT, by clicking on it. Set both the Range and the 
Display Range to the vector [2 15]. 

> Similarly name the triangles as low, medium and high with parameters [4 5], [4 5 
6] and [5 6] 

> Select the output variable, GT, by clicking on it. Set both the Range and the 
Display Range to the vector [90 300]. 

> Select Add MFs. from the Edit menu and add five trimf curves to the output 
variable GT. 

> Similarly name the triangles as very low, medium high and very high with 
parameters [140 160], [140 160 180], [160 180 200], [180 200 220] and [200 220] 

Now the final system should look something like the fig 7.7 below, 
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Fib' Variables 
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Type output 
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Current Membership Function (click on MF to select) 
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Display Range j [90 300] 

Help Close 


Selected variable "GT” 


Fig 7.7 Final FIS System for Grain Temperature (GT) 


7.3.4 Rule Editor 

Now one is ready to write down the rules as give in the table 7.2. To call up the 
Rule Editor, go to the View menu and select Edit rules, or type ruleedit at the command 
line. 

Constructing rules using the graphical Rule Editor interface as shown in fig 7.8 
below is fairly self-evident. Based on the descriptions of the input and output variables 
defined with the FIS Editor, the Rule Editor allows you to construct the rule statements 
automatically, by clicking on and selecting one item in each input variable box, one item 
in each output box, and one connection item. Choosing none as one of the variable 
qualities will exclude that variable from a given rule. Choosing not under any variable 
name will negate the associated quality. Rules may be changed, deleted, or added, by 
clicking on the appropriate button. 
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Fig 7.8 Rule Editor 


The Rule Editor also has some familiar landmarks, similar to those in the FIS 
Editor and the Membership Function Editor, including the menu bar and the status line. 
The Format pop-up menu is available from the Options pull-down menu from the top 
menu bar; this is used to set the format for the display. Now one has to stait entering the 
rules as per the table 7.2 by selecting the entries from the boxes. Once all the 27 rules 
have been entered the GT.fis system is fully ready for viewing the rules and surfaces, 
which can be done using the view menu. 

7.3.5 Rule Viewer and Surface Viewer. 

On selecting the rule viewer the fig7.9 pops up as under, 
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4 Rule Viewer: GT1 



Fig 7.9 Rule Viewer 


The Rule Viewer allows you to interpret the entire fuzzy inference process at 
once. The Rule Viewer also shows how the shape of certain membership functions 
influences the overall result. Since it plots every part of every rule, it can become 
unwieldy for particularly large systems, but, for a relatively small number of inputs and 
outputs, it performs well (depending on how much screen space you devote to it) with up 
to 30 rules and as many as 6 or 7 variables. 

The Rule Viewer shows one calculation at a time and in great detail. In this sense, 
it presents a sort of micro view of the fuzzy inference system. If you want to see the 
entire output surface of your system, that is, the entire span of the output set based on the 
entire span of the input set, you need to open up the Surface Viewer. This is the last of 
the five basic GUI tools in the Fuzzy Logic Toolbox, and you open it by selecting View 
surface... from the View menu. The following screen pops up as in fig 7.10 below. 
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Fig 7.10 Surface Viewer Plot for CTI & GHDT Vs GT 

Upon opening the Surface Viewer, we are presented with three-dimensional plots 
that MATLAB can adeptly manage. When we move beyond three dimensions as in this 
case overall, we start to encounter trouble displaying the results. Accordingly, the Surface 
Viewer is equipped with pop-up menus that let you select any two inputs and any one 
output for plotting. Just below the pop-up menus are two text input fields that let you 


i Surface Viewer: GT1 


HU El 


determine how many x-axis and y-axis grid lines you want to include. This allows you to 
keep the calculation time reasonable for complex problems. Pushing the Evaluate button 
initiates the calculation, and the plot comes up soon after the calculation is complete. To 
change the x-axis or y-axis grid after the surface is in view, simply change the appropriate 
text field, and click on either X-grids or Y-grids, according to which text field you 
changed, to redraw the plot. 
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7.4 Final Output and Results 

Similarly the remaining sub systems were developed as MIX_TEMP, PMT, GBD, 
COKEDPROCITY, CCS, HMOR, SCR and OXYDATION RESISTANCE. The membership 
functions for each of these are attached as per Fig 7.1 1 to 7.18 of Appendix G. The rule sets for 
each of these are attached as Table 4.3 to Table 4.10 in Appendix A. The surface plots up to the 
Green BD stage are attached as per Fig 7.19 to 7.39 of Appendix H. The first 4 sub systems were 
combined using a matlab programme file called GBD_LONG.m. This programme accepts input 
from files and the out of the first is fed to the next along with additional inputs and this chain 
goes on till we get the final out as Green Bulk Density. Individual matlab programme files have 
also been made to accept inputs in vector form. In addition programmes have also been made so 
that the user can input the data file from the keyboard when prompted to do so. The results 
obtained using sample data are attached as per tables 7.1 1 to 7.15 of Appendix I. 

Validation of the expert based system requires a large volume of low noise data so that 
the membership functions can be modified to achieve the desired results. The non availability 
such a data set has meant that this has to be vigorously followed up in the future and the Expert 
based Fuzzy Inference system should be fine tuned to achieve best results. 
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CHAPTER 8 


DEVELOPMENT OF FIS FORM PLANT DATA 
8.1 Overview of ANFIS SYSTEM 

8.1.1Sugeno-Type Fuzzy Inference 

The fuzzy inference process we’ve been referring to so far is known as 
Mamdani’s fuzzy inference method. It’s the most commonly seen fuzzy methodology. In 
this section we discuss the Sugeno method of fuzzy inference as the ANFIS EDITOR of 
matlab toolbox utilizes this method only. The first two parts of the fuzzy inference 
process, fuzzifying the inputs and applying the fuzzy operator, are exactly the same. The 
main difference between Mamdani-type of fuzzy inference and Sugeno-type is that the 
output membership functions are only linear or constant for Sugeno-type fuzzy inference. 

A typical fuzzy rule in a zero-order Sugeno fuzzy model which has the form if x is 
A and y is B then z = k, where A and B are fuzzy sets in the antecedent, while A: is a crisply 
defined constant in the consequent. When the output of each rule is a constant like this, 
the similarity with Mamdani’s method is striking. The only distinctions are the fact that 
all output membership functions are singleton spikes, and the implication and aggregation 
methods are fixed and can not be edited. The implication method is simply 
multiplication, and the aggregation operator just includes all of the singletons. 

* Advantages of the Sugeno Method 

• It’s computationally efficient. 

• It works well with linear techniques (e.g., PID control). 

• It works well with optimization and adaptive techniques. 

• It has guaranteed continuity of the output surface. 

• It’s well-suited to mathematical analysis. 

* Advantages of the Mamdani Method 

• It’s intuitive. 

• It has widespread acceptance. 

• It’s well-suited to human input. 

U Anfis and the ANFIS Editor GUI 

In this section we discuss the use of the function anfis and the ANFIS Editor GUI in the 
? uzzy Logic Toolbox. These tools apply fuzzy inference techniques to data modeling. As you 
ave seen from the other fuzzy inference GUIs, the shape of the membership functions depends 
n parameters, and changing these parameters will change the shape of the membership function. 
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Instead of just looking at the data to choose the membership function parameters, we will see 
how membership function parameters can be chosen automatically using these Fuzzy Logic 
Toolbox applications. 

8.2.1 Data Modeling 

For modeling the manufacturing process of RSP we were in the initial stages 
given a data set of 86 observations consisting of 10 inputs and 1 output, the same is in 
table 4.13 attached as Appendix C. An additional set of 31 data sets (in table 4.12 of 
Appendix C) was sent subsequently making a total of 117 observations hence the 
modeling was done for both the 86 and 117 data sets separately. Since the membership 
functions chosen for the various attributes in the expert system developed by expert 
advice was to certain extent arbitrary. Rather than choosing the parameters associated 
with a given membership function arbitrarily, these parameters could be chosen so as to 
tailor the membership functions to the input/output data in order to account for these 
types of variations in the data values. This is where the neuro-adaptive learning 
techniques incorporated into anfis in the Fuzzy Logic Toolbox can help. The inputs and 
output for RSP for which data was made available is as under, 

♦ Inputs 

> Combustion chamber temperature 

> Grain heater drum temperature 

> Grain retention time 

> Mixer temperature 

> Mixing time 

> Mix temperature 

> Storage time 

> Rolls 

> Press mix temperature 

> Forming pressure 

♦ Outputs 

> Green bulk density (GBD) 

8.2.2 Learning and Inference through ANFIS 
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The basic idea behind these neuro-adaptive learning techniques is that these 
techniques provide a method for the fuzzy modeling procedure to learn information about 
a data set, in order to compute the membership function parameters that best allow the 
associated fuzzy inference system to track the given input/output data. This learning 
method works similarly to that of neural networks. The Fuzzy Logic Toolbox function 
that accomplishes this membership function parameter adjustment is called anfts 
( adaptive neuro-fuzzy inference System). Using a given input/output data set, the toolbox 
function anfis constructs a fuzzy inference system (FIS) whose membership function 
parameters are tuned (adjusted) using either a back propagation algorithm alone, or in 
combination with a least squares type of method. This allows your fuzzy systems to learn 
from the data they are modeling. 

8.2.3 Importance of Accurate Data for Good Modeling 

The modeling approach used by anfis is similar to many system identification 
techniques. First, you hypothesize a parameterized model structure (relating inputs to 
membership functions to rules to outputs to membership functions, and so on). Next, you 
collect input/output data in a form that will be usable by anfis for training. You can then 
use anfis to train the FIS model to emulate the training data presented to it by modifying 
the membership function parameters according to a chosen error criterion. In general, this 
type of modeling works well if the training data presented to anfis for training 
(estimating) membership function parameters is fully representative of the features of the 
data that the trained FIS is intended to model. This is not always the case as in the data 
given to us by RSP has only 117 observations where we would require about 3000 to 
4000 data sets. More over the data collected had noisy measurements, and hence the 
training data could be representative of all the features of the data that will be presented 
to the model. 

8.2.4 Model Validation Using Testing Data and Checking Data Sets 

Model validation is the process by which the input vectors from input/output data 
sets on which the FIS was not trained, are presented to the trained FIS model, to see how 
well the FIS model predicts the corresponding data set output values. This is 
accomplished with the ANFIS Editor GUI using the testing data set. In this case the 31 
data set received subsequently was used as testing data set. The basic idea behind using a 
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checking data set for model validation is that after a certain point in the training, the 
model begins over fitting the training data set. In principle, the model error for the 
checking data set tends to decrease as the training takes place up to the point that over 
fitting begins, and then the model error for the checking data suddenly increases. 

8.2.5 The ANFIS Editor GUI 

To get started with the ANFIS Editor GUI, type anfisedit on the command 
window and the following window as in fig8.1 will pop up. 


8.2.5.1 Loading Data 

From the GUI above one can, Load data (training and testing, and checking) by 
selecting appropriate radio buttons in the Load data portion of the GUI and then 
selecting Load Data the loaded data is plotted on the plot region. Here we load the 86 
data sets for training and 3 1 data sets for testing. Since the numbers of data sets were not 
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FIG 8.1 Anfis GUI 



sufficient checking data could not be loaded. This is how it looks once the data has been 


loaded as in fig 8.2 
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Fig 8.2 Loading of Training and Testing Data 
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8.2.S.2 Generating the FIS Model 

Geneiate an initial FIS model or load an initial FIS model using the 
options in the Generate FIS portion of the GUI. 

y Flere we generate two Fis models using Grid partitioning option where we 
specify 3 MFs for each input and 3 MFs for the outputs. The two Fis is for 
the data sets broken up into two parts the first having 5 input and 1 out put 
and the second having the output of the first as one of its input along with 
remaining 4 of the 10 inputs. This splitting of the inputs had to be done 
because with 10 inputs taken together and each with 3 MFs, the system 
would have to generate 3 10 rules and hence was running out of memory. 
Further the membership functions selected were triangular and trapezoidal 
only, because for all other memberships’ functions (gbell, Gaussian, 
sigma, dsigma, psigma, pi, S and Z membership functions) the results 
were not very encouraging. Hence two sets of Fis were generated, one set 
for triangular membership functions and the other for trapezoidal. 


Anfis Model Structure Hill El 



Fig 8.3 ANFIS Structure 
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K View the FIS model structure once an initial FIS has been generated or 
loaded by selecting the Structure button. This is how the structure looks 
like as in fig 8.3 above 
8.2.5.3 Training the Generated FIS 

> Train the Fis model using the parameter optimization method as Hybrid 
(mixture of back propagation and least squares) and number of epochs as 
300 and margin of error as 0.01 and click train now. As the training 
progresses the following screen is displayed as per fig 8.4 below. And 
once training is completed the screen displayed is shown in fig 8.5. 
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Fig 8.4 Training of Data for Generation of FIS 





4 Anfis Editor. mixtemp_86_3MF_TRIANGLE 


Generate FIS 


Train FIS 


Optim. Method: 
[hybrid 

Error Tolerance: 


Plot against: 

<* Training data 
C Testing data 


<* Training 


^ Load from worksp. 
(* Grid partition 
C Sub. clustering 


r Testing 


C worksp. 


C Checking 


r Checking data 


Generate FIS 


Fig 8.5 Training Completed and Training Error Saturates At 1.2582 


> Once the Fis have been generated save them to the disk, as has been done 
in the case of RSP data of 86 and 117 data sets. A total of 8 Fis have been 
developed as per table 8.3 below. 

8.2.6 Viewing/Editing the Mix Temperature Fis System Generated By ANFIS 

Once the Fis system have been saved to the disk type fuzzy on the command 
window to open the fuzzy editor, and from the file menu open the Mix_Temp Fis and you 
will see the structure as in fig 8.6. In this one has to rename the input and output 
parameters as desired as by default they are named as input 1, input2 and so on, in our 
case they have been renamed as CTI, GHDT, RT, MIXERTEMP, MIXINGTIME and 
MIXTEMP. The procedure of naming the parameters has already been explained in 
chapter 7. 
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NAME OF 

FIS 

NO OF 

DATA 

POINTS 

TYPE OF 

MEMBERSHIP 

FUNCTION 

INPUTS WITH 

3 MF EACH 

OUTPUT WITH 

3 MF 

MIXJTEMP 

86/117 

TRIANGULAR 

CTI 

MIX_ TEMP 




GHDT 





RT 





MIXER TEMP 





MIXING TIME 


GBD 

86/117 

TRIANGULAR 

MIX_ TEMP 

GBD 




ST 





ROLLS 





PMT 





FP 


MIXJTEMP 

86/117 

TRAPEZOIDAL 

CTI 

MIX_ TEMP 




GHDT 





RT 





MIXER TEMP 





MIXING TIME 


GBD 

86/117 

TRAPEZOIDAL 

MIX_ TEMP 

GBD 




ST 





ROLLS 





PMT 





FP 



Table 8.3 Total of 8 Fis Systems Developed For RSP 
8.2.6. 1 Viewing/Editing the Membership Functions Generated By ANFIS 

In order to view/edit the membership functions one selects edit 
membership functions and the following window will pop up showing all the 
membership functions that were generated for the input and output variables. As 


84 




Fig 8.6 Fis System for Mix Temp 

a sample only the membership functions for CTI has been shown in fig 8.7 below. 
The membership functions for the remaining parameters GHDT. Retention Time, 
Mixer Temperature, Mixing time, are shown as figures 8.14 to 8.17 per Appendix 
K 

8.2.6.2 Viewing/Editing the Rules Generated By ANFIS 

For this one has to select edit rules and the following window is displayed 
as shown in fig 8.8 below. A total of 3 5 = 243 rules are generated by the ANFIS 
editor. In case any of the rules are to be changed or modified they can be done as 
per procedure explained in chapter 7. In the case of RSP data model no rules have 
been modified as the results obtained were very encouraging. It may be noted that 
just by changing the variable names and the membership function names the rules 
get automatically updated to the new names. 
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Fig 8.7 Membership functions for CTI generated by ANFIS editor 
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Fig 8.8 Rules Generated By ANFIS Editor Total Of 3 5 = 243 Rules 
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8.2.6.3 Rule Viewer and Surface Viewer 


The rule viewer and the surface viewer as already brought are read only 
windows and no editing can be done using them. But they give a bird’s eye view 
of the relation between the inputs and out put. The windows generated on 
selecting these options are shown in figures 8.9 and 8.10. For the surface viewer 
only one of the combinations i.e. CTI and GHDT plotted against Mix 
temperature has been shown in the chapter the remaining have been attached as 
figures 8.18 to 8.20 of Appendix L. 

8.2.7 Viewing/Editing the GBD Fis System Generated By ANFIS 

On similar lines the Fis system for GBD has be generated by following the above 
steps and the results and out puts obtained have been attached as figures 8.21 to 8.32 in 
Appendix M. 


4 Rule Viewer. miKtemp_86_3MF_TRlANGLE 


File Edit View Options 


GHDT ■ 504 RT * 6.5 MIXERTEMP = 113 MIXINGTIME = 18.5 MIXTEMP = 202 



Fig 8.9 Rule Viewer for Mix temperature 
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Fig 8.10 Surface Viewer for CTI and GHDT Plotted Against Mix Temperature 

8.3 Validation Of Training Data Against Testing Data For 86 Data Set Vs 31 Data Set. 

Once the training and testing data sets, Fis have been loaded from the disk now 
the test Fis radio buttons can be used to test the system for which Plot against testing data radio 
button is selected. The result obtained is shown in the figure 8.11 below. As can be seen from the 
plot which shows the desired out put in blue dots and the actual out put obtained from the Fis 
system trained by 86 data set with red crosses, that the out put obtained using the testing data has 
far to high an average error of 78.1146, which confirms our contention that with just 86 data sets 
the system will not be able to capture the entire feature space of the process for which we require 
at least 3000 to 4000 data sets. 

8.4 Results Obtained Using 117 Data Sets for Mix temperature Fis 

The Fis systems generated using the 117 data set was tested and checked with the 86 and 
31 data sets respectively. For this the test Fis was plotted against testing data and checking data 
and the results are shown in figure 8.12 and 8.13 respectively. It can be seen that the average 
error obtained has drastically dropped to 2.0939 and 2.9476 respectively. Which further 
substantiates the importance of training with larger volume of data? 
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Fig 8.11 Testing the Fis System against Testing 
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Fig 8.12 Results using Fis trained withll7 data sets Plot against tes g 
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Fig 8.13 Results using Fis trained with!17 data sets Plot against checking data of 31 sets 
8.5 The Final Results for the Combined Systems 

The Fis generated for Mix temperature and GBD were combined as explained earlier in 
this chapter using a matlab programme for the same. The programmes were also made to accept 
input for the user in vectors from or as a complete file from the keyboard itself. The outputs 
obtained for 86 data set Fis and That for 117 data set Fis are shown in table 8.4 and 8.5 attached 
as Appendix J. The 86 data set system results for Triangular membership and Trapezoidal 
membership have been compared with the Expert system developed in prolog as well. The out 
put of the systems developed for 117 data sets has been compared for Triangular membership 
and Trapezoidal membership. It has been observed that the Triangular system has more accurate 
and stable out put compared to that of trapezoidal system in both cases. When evei further new 
data sets get collected they should be test and checked with the 117 data set system for 
Validation as explained vide paragraph 8.3 and in case the results are satisfactory the system 
need not be trained afresh, but as and when the variation in data collected with minimum noise 
increases it will be better to retrain the system so that it is able to fully cover the feature space of 
the data distribution. 
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CHAPTER 9 


COMPARISON BETWEEN THE MATLAB AND PROLOG RASED 

FUZZY INFERENCE SYSTEMS 

9.1 Prolog based system 

The prolog based FIS system is not equipped with the GUI and hence is not very 
user friendly for the plant operators [10]. It requires greater training for the operator to 
use the system effectively. 

The prolog based system is not equipped with the surface view and rule viewer 
graphics that give an intuitive feel of the behavior of the input parameters and the rules 
that fire the system. This has its advantages as well in that the prolog system has a faster 
response to the input data as compared to the inputs to the Matlab system; this can have 
an effect in the future as the number of data collected increases many fold with the 
passage of time. 

The prolog based system has only three types of membership functions, including 
the user defined function, compared to the matlab based system which can cater for 
eleven types of membership functions, excluding the user defined membership functions, 
as already brought out in chapter 7. However in the application for RSP only two types of 
the memberships could be used. 

The implication method for prolog is only Mamdani, where as matlab has both 
Mamdani and Sugeno type and in addition user defined implications can also be used. 

The Anfis editor of matlab is unique in that it generates the membership functions 
and rules on its own to fit the data input and has a built in safe guard to avoid over fitting. 
The only condition is that the data should be of low noise and of sufficient quantity to be 
able to be broken up into, separate training, teat and checking data. No such provision 
exists in the prolog system. All rule and membership functions have to be generated by 
the user/programmer himself. 

The coding of programmes for matlab is very simple and user friendly and any 
changes to the programme can be implemented with ease and minimum amount of fuss. 
More over debugging of the programmes are very simple with extensive literature and 
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online help icadily available lor any needed assistance/help. The same can not be said of 
the prolog software. 

The results obtained on test data has shown that the performance of both systems 
is satisfactory and hence result wise both are equally good. 

9.2 Future Recommendation’s 

The Validation of both the systems developed has been held up for want of low 
noise data in sufficient quantity. This data should be made available as early as possible 
so as to fine tune the system either by modifying the membership functions or the rules 
base or both as the case may be. This may not be a very major job but is a very important 
one as only then can the system be put to any effective use. 

The backward chaining systems developed in the Bayesian and FIS engines used 
for fault diagnoses have also to fine tuned for want of proper rules sets as they have not 
yet been provided by the experts nor is any tangible data available on the matter to 
generate these rules using tools such as Cart, ID3 etc. 

The brick manufacturing process at RSP has a provision for it to be run on 
automation, which has not been done so far. The FIS system developed should be 
integrated to the automated system whenever it is made functional, only the will it be 
really possible to achieve the final target Brick life that the factory has in mind. 

A graphic user interface should be developed for both the matlab and prolog 
systems so that the operator is able to run the desired modules in the programme by just a 
click of a button from his desktop. 

9.3 Conclusion 

In the mean while the FIS system developed can be used to train the operators on 
its use, so as to understand its capability and limitations. The much needed data to fine 
tune the system should be collected as early as possible for further processing. The 
control of the input parameters can still be achieved in the limited sense that the system 
has been able to learn form the however limited data that has already been used to 
develop and train the FIS system. The Expert based FIS can be fully implemented at this 
stage to monitor its performance even beyond the GDB stage. 
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Appendix A 


EXPERT DASED DATA SET FOR RULE GENERATION 
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Table 4.2 for Grain Temperature & t he ranges for High, Medium & Low 
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Table 4.3 for Mix Temperature & the ranges for Hifrh, Medium & Low 
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VERY HIGH 

VERY LOW 

LOW 

MEDIUM 

MEDIUM 

HIGH 

VERY LOW 

LOW 

MEDIUM 

MEDIUM 

MEDIUM 

LOW 

LOW 

MEDIUM 

MEDIUM 

LOW 

MEDIUM 

LOW 

MEDIUM 

MEDIUM 

VERY LOW 

. MEDIUM; .: 

LOW 

MEDIUM 

LOW 

VERY HIGH 

VERY LOW . 

LOW 

MEDIUM 

LOW 

HIGH 


LOW 

MEDIUM 

LOW 

MEDIUM 

i' \ 1 J ; 'Iff 


MEDIUM 

LOW 

LOW 

", „ ' '1 

LOW 

MEDIUM 

LOW 

VERY LOW 

' ' ’ ■ 11 ‘ ; v 

LOW 

LOW 

HIGH 

VERY HIGH 

ilPIfUS 

LOW 

LOW 

HIGH 

HIGH 


LOW 

LOW 

HIGH 

MEDIUM 


, LOW 

LOW 

HIGH 

LOW 


LOW 

LOW 

HIGH 

VERY LOW 


■EM 


MEDIUM 


* u , i - ' . 1111 » 

IBlifli 

LOW 

MlaMUJSl 



’' r m 

LOW 

MEDIUM 

MEDIUM 

;I~ pi 


LOW 

MEDIUM 

LOW 

■ 

fei 

IliS 

, 

■ . 

-im 

LOW 

MEDIUM 

VERY LOW 



LOW 

LOW 

VERY HIGH 



LOW 

LOW 

HIGH 

, v ■; . Lu Mpl 

LOW 

LOW 

LOW 

MEDIUM 



LOW 

LOW 

LOW 


LOW 

LOW 

LOW 

LOW 

v is < v 


TABLE4.5 for GBP & the ranees for Hiflh. Medium & 


LOW 


HIGH 


HIGH 


HIGH 


LOW 


HIGH 


HIGH 


MEDIUM 


ow{2gf 2) 
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f PITCH 

T TEMI’ERIN 

CCS 

t high 

J PERFECT 

VERYHiGlT 

f HIGH 

j UNDER 

"veryjrgfT 

1 HIGH 

| TT 

~1gh ' 

(MEDIUM 

j perfect 

W®gh“ 

f MEDIUM 

] Under" 

TnoT" 

f MEDIUM 

OVER 

Tigh 

^ LOW 

PERFECT 

HIGH 

! LOW 

“under 

MEDIUM 

j LOW 

OVER 

VERYJHIGH 

HIGH 

PERFECT^ 

VERYJHIGH 

HIGH 

UNDER 

HIGH 

f HIGH OVER 

HIGH 

(MEDIUM 

PERFECT^ 

HIGH 

MEDIUM 

undeT 

MEDIUM 

"MEDIUM j 

OVER 

MEDIUM 

f LOW 

PERFECT^ 

MEDIUM 

LOW 

UNDER 

MEDIUM 

LOW | 

OVER 

MEDIUM 

HIGH j PERFECT^ 

HIGH 

HIGH J UNDER 

HIGH 

HIGH j OVeT^ 

MEDIUM 

MEDIUM 

' MEDIUM j PERFECT 

“medium! undeT^ 

MEDIUM 

medium] oveT^ 

LOW 

j LOW PERFECT 

LOW 

LOW UNDER 

MEDIUM 

r low 

0VER 

LOW 

HIGH 

PERFECT 

MEDIUM 

HIGH 

UNDER 

MEDIUM 

HIGH 

j OVER __ 

LOW 

Tedium 

(“medium 

PERFECT 

pMEDIUM - 

UNDER ____ 

"“medium 

“MEINUmI 

OVER^T 

LOW 

low perfecTI 

LOW 

LOW ( 

UNDER ____ 

LOW 

VERY_LOW 

'Tow | 

OVER 

HIGH T 

PERFECT 

VERY_LOW 


JL JUVJ A A ~ — — 

TARTF4.6 for Rules'oS^^ 


98 


PITCH 

TEMPERING 

CCS 

HIGH 

UNDER 

LOW 

HIGH 

OVER 

VERYJLOW 

MEDIUM 

PERFECT 

— _ 

VERY_LOW 

MEDIUM 

UNDER 

LOW 

[medium 

OVER 

VERYJLOW 

LOW 

PERFECT 

VERY_LOW 

LOW 

UNDER 

VERY_LOW 

LOW 

OVER 

VERYJLOW 


TABLE4.6 for Rules of CCS(2 of 2) 


GREEN BD 

PITCH 

QUANTITY 

COKED 

POROSITY 

VERY HIGH 

HIGH 

MEDIUM 

VERY HIGH 

MEDIUM 

LOW 

VERY HIGH 

LOW 

LOW 

HIGH 

HIGH 

HIGH 

HIGH 

MEDIUM 

MEDIUM 

HIGH 

LOW 

LOW 

MEDIUM 

HIGH 

HIGH 

MEDIUM 

MEDIUM 

MEDIUM 

MEDIUM 

LOW 

MEDIUM 

LOW 

HIGH 

HIGH 

LOW 

MEDIUM 

MEDIUM 

LOW 

LOW 

MEDIUM 

VERY LOW 

HIGH 

HIGH 

VERY LOW 

MEDIUM 

HIGH 

VERYJ.OW 

LOW 

HIGH 


TABLE4.7 for Rules of CP 


CARBON 

CONTENT 

T ~ — ' J 

] 

j 

GREEN BD 

METAL 

POWDER 

COLD 

CRUSHING 

STRENGTH 

HMOR 

VERY HIGH 

LOW 

4- — - — 

HIGH 

HIGH 

VERY VERY HIGH 

VERY HIGH 

LOW 

HIGH 

MEDIUM 

VERY VERY HIGH 

VERY HIGH 

LOW 

t— — — — 

HIGH 

LOW 

VERY VERY HIGH 

VERY HIGH 

LOW 

. . . ... ... 

MEDIUM 

HIGH 

VERY VERY HIGH 

VERY HIGH 

LOW 

j — — 

MEDIUM 

MEDIUM 

VERY VERY HIGH 

VERY HIGH 

LOW 

MEDIUM 

LOW 

VERY HIGH 

VERY HIGH 

LOW 

LOW 

HIGH 

HIGH 

VERY HIGH 

LOW 

LOW 

MEDIUM 

HIGH 

VERY HIGH 

LOW 

LOW 

LOW 

MEDIUM 

VERY HIGH 

VERY.. LOW 

HIGH 

HIGH 

VERY VERY HIGH 

VERY HIGH 

VERY LOW 

HIGH 

MEDIUM 

VERY HIGH 

VERY HIGH 

VERY LOW 

HIGH 

LOW 

VERY HIGH 

VERY HIGH 

VERYLOW 

MEDIUM 

HIGH 

HIGH 

VERY HIGH 

VERY_LOW 

MEDIUM 

MEDIUM 

HIGH 

VERY HIGH 

VERY LOW 

MEDIUM 

LOW 

HIGH 

VERY HIGH 

VERY LOW 

LOW 

HIGH 

MEDIUM 

VERY HIGH 

VERYJ.OW 

LOW 

MEDIUM 

MEDIUM 

VERY HIGH 

veryTow. 

LOW 

LOW 

MEDIUM 

HIGH 

MEDIUM 

HIGH 

HIGH 

VERY_ VERY HIGH 

HIGH 

MEDIUM 

HIGH 

MEDIUM 

VERY VERY ..HIGH 

HIGH 

MEDIUM 

HIGH 

LOW 

VERY VERY HIGH 

HIGH 

MEDIUM 

MEDIUM 

HIGH 

VERY VERY HIGH 

HIGH 

MEDIUM 

MEDIUM 

MEDIUM 

VERY VERY HIGH 

HIGH 

MEDIUM 

MEDIUM 

LOW 

VERY HIGH 

HIGH 

MEDIUM 

LOW 

HIGH 

HIGH 

HIGH 

MEDIUM 

LOW 

MEDIUM 

MEDIUM 

HIGH 

MEDIUM 

LOW 

LOW 

MEDIUM 

HIGH 

LOW 

HIGH 

HIGH 

VERY VERY HIGH 

HIGH 

LOW 

HIGH 

MEDIUM 

VERY VERY HIGH 

HIGH 

LOW 

HIGH 

LOW 

VERY HIGH 

HIGH 

LOW 

MEDIUM 

HIGH 

VERY HIGH 

HIGH 

LOW 

MEDIUM 

MEDIUM 

VERY HIGH 

HIGH 

LOW 

MEDIUM 

LOW 

HIGH 

HIGH 

LOW 

LOW 

HIGH 

HIGH 


TABLE4.8 RULES FOR HMORIl of 31 


100 



HIGH 

LOW 

LOW 

MEDIUM 

MEDIUM 

HIGH 

I OW 

LOW 

LOW 

MEDIUM 

HIGH 

VERYJLOW 

HIGH 

HIGH 

VERY HIGH 

HIGH 

MERYLOW 

HIGH 

MEDIUM 

VERY HIGH 

HIGH VKRYJ.OW 

HIGH 

LOW 

VERY HIGH 

HIGH V ; H R Y LOW 

MEDIUM 

HIGH 

HIGH 

HIGH VHRYLOW 

MEDIUM 

MEDIUM 

HIGH 

HIGH V T HKY LOW 

MEDIUM 

LOW 

MEDIUM 

HIGH WHRY LOW 

LOW 

HIGH 

MEDIUM 

HIGH [ VERY LOW 

LOW 

MEDIUM 

MEDIUM 

HIGH J VERY LOW 

LOW 

LOW 

MEDIUM 

WED IUIM f MEDIUM 

HIGH 

HIGH 

VERY VERY HIGH 

NED HUM 

MEDIUM 

HIGH 

MEDIUM 

VERY VERY HIGH 

NEDIUN 

MEDIUM 

HIGH 

LOW 

VERY VERY HIGH 

IMEDIUN 

MEDIUM 

MEDIUM 

HIGH 

VERY VERY HIGH 

NEDIUN 

MEDIUM 

MEDIUM 

MEDIUM 

VERY HIGH 

NEDIUN 

MEDIUM 

MEDIUM 

LOW 

HIGH 

NEDIUN 

MEDIUM 

LOW 

HIGH 

MEDIUM 

NEDIUN 

MEDIUM 

LOW 

MEDIUM 

MEDIUM 

NEDI UN 

MEDIUM 

LOW 

LOW 

LOW 

NEDI UN 

LOW 

HIGH 

HIGH 

VERY HIGH 

MEDI UN 

LOW 

HIGH 

MEDIUM 

VERY HIGH 

MEDI UN 

LOW 

HIGH 

LOW 

VERY HIGH 

MEDI UN 

LOW 

MEDIUM 

HIGH 

VERY HIGH 

MEDI DM 

LOW 

MEDIUM 

MEDIUM 

VERY HIGH 

MEDIUM 

LOW 

MEDIUM 

LOW 

HIGH 

MEDIUM 

LOW 

LOW 

HIGH 

MEDIUM 

MEDIUM 

LOW 

LOW 

MEDIUM 

MEDIUM 

MEDIUM 

LOW 

LOW 

LOW 

LOW 

MEDIUM 

VERY LOW 

HIGH 

HIGH 

HIGH 

MEDIUM 

VERY LOW 

HIGH 

MEDIUM 

MEDIUM 

MEDIUM 

VERY LOW 

HIGH 

LOW 

MEDIUM 

MEDIUM 

VERY LOW 

MEDIUM 

HIGH 

MEDIUM 

MEDIUM 

VERY LOW 

MEDIUM 

MEDIUM 

MEDIUM 

_ MEDIUM 

VERY LOW 

MEDIUM 

LOW 

LOW 

MEDIUM 

VERY LOW 

LOW 

HIGH 

LOW 

medium 

YERY LOW 

LOW 

MEDIUM 

VERY LOW 
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r 

GREEN BD 

CARBON 

CONTENT 

METAL 

POWDER 

OXIDATION 

resistance 

VERYHIGH 

VERY LOW 

HIGH 

MEDIUM 

VERY HIGH 

VERY LOW 

MEDIUM 

LOW 

; VERY HIGH 

VERY LOW 

LOW 

VERY LOW 

HIGH 

VERY LOW 

HIGH 

LOW 

HIGH 

VERY LOW 

MEDIUM 

VERY LOW 

HIGH 

VERY LOW 

LOW 

VERY LOW 

HIGH 

LOW 

HIGH 

MEDIUM 

HIGH 

LOW 

MEDIUM 

LOW 

HIGH 

LOW 

LOW 

LOW 

MEDIUM 

HIGH 

HIGH 

VERY HIGH 

MEDIUM 

HIGH 

MEDIUM 

HIGH 

MEDIUM 

HIGH 

LOW 

MEDIUM 

MEDIUM 

MEDIUM 

HIGH 

HIGH 

MEDIUM 

MEDIUM 

MEDIUM 

HIGH 

_ MEDIUM 

MEDIUM 

LOW 

MEDIUM 

MEDIUM 

LOW 

HIGH 

MEDIUM 

MEDIUM 

LOW 

MEDIUM 

LOW 

MEDIUM 

LOW 

^ Tir . r T - ^.- nr r _ ... 

LOW 

LOW 

MEDIUM 

VERY LOW 

HIGH 

LOW 

MEDIUM 

VERY LOW 

MEDIUM 

LOW 

MEDIUM 

VERY LOW 

LOW 

VERY_LOW 

LOW 

VERY HIGH 

HIGH 

VERY HIGH 

LOW 

VERY HIGH 

MEDIUM 

VERY HIGH 

LOW 

VERY HIGH 

LOW 

HIGH 

LOW 

VERY LOW 

HIGH 

LOW 

LOW 

VERY LOW 

MEDIUM 

VERYJLOW 

LOW 

VERY LOW 

LOW 

VERY LOW 

LOW 

HIGH 

HIGH 

HIGH 

LOW 

HIGH 

MEDIUM 

HIGH 

LOW 

HIGH 

LOW 

MEDIUM 

LOW 

MEDIUM 

HIGH 

MEDIUM 

LOW 

MEDIUM 

MEDIUM 

LOW 
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LOW 

MEDIUM 

LOW 

LOW 

LOW 

LOW 

HIGH 

LOW 

LOW 

LOW 

MEDIUM 

VERY LOW 

LOW 

LOW 

LOW 

VERY LOW 

VERY LOW 

HIGH 

HIGH 

MEDIUM 

VERY LOW 

HIGH 

MEDIUM 

MEDIUM 

VERYJLOW 

HIGH 

LOW 

LOW 

VERY. LOW 

MEDIUM 

HIGH 

MEDIUM 

VERY_LOW 

MEDIUM 

MEDIUM 

LOW 

VERY., LOW 

MEDIUM 

LOW 

LOW 

VERYJLOW 

VERY. HIGH 

HIGH 

VERY HIGH 

VERY LOW 

VERY HIGH 

MEDIUM 

HIGH 

VERY LOW 

VERY HIGH 

LOW 

MEDIUM 


TABLE4.9 RULES FOR OXIDATION RESISTANCE (2 of 21 


CARBON 

CHEMICAL 

COKED 

SLAG CORROSION 

CONTENT 

PURITY 

POROSITY 

RESISTANCE 

MEDIUM 

HIGH 

HIGH 

MEDIUM 

MEDIUM 

HIGH 

MEDIUM 

HIGH 

MEDIUM 

HIGH 

LOW 

HIGH 

MEDIUM 

MEDIUM 

HIGH 

MEDIUM 

MEDIUM 

MEDIUM 

MEDIUM 

MEDIUM 

MEDIUM 

MEDIUM 

LOW 

HIGH 

MEDIUM 

LOW 

HIGH 

LOW 

MEDIUM 

LOW 

MEDIUM 

LOW 

MEDIUM 

LOW 

LOW 

LOW 

LOW 

HIGH 

HIGH 

LOW 

LOW 

HIGH 

MEDIUM 

MEDIUM 

LOW 

HIGH 

LOW 

MEDIUM 

LOW 

MEDIUM 

HIGH 

LOW 

LOW 

MEDIUM 

MEDIUM 

MEDIUM 

LOW 

MEDIUM 

LOW 

MEDIUM 

LOW 

LOW 

HIGH 

VERY_LOW 

LOW 

LOW 

MEDIUM 

VERYLOW 

LOW 

LOW 

LOW 

LOW 
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CARBC )N 
CONTENT 

( HLMICAL 
PURITY 

COKED 

POROSITY 

SLAG CORROSION 
RESISTANCE 

VERY i ( )W 

HIGH 

HIGH 

LOW 

VERY I <)W 

HIGH 

MEDIUM 

MEDIUM 

VERY LOW 

HIGH 

LOW 

MEDIUM 

VERY LOW 

MEDIUM 

HIGH 

VERY_LOW 

VERY LOW 

MEDIUM 

MEDIUM 

LOW 

VERY LOW 

MEDIUM 

LOW 

LOW 

VERY LOW 

LOW 

HIGH 

VERYJLOW 

VERY LOW 

LOW 

MEDIUM 

VERYJLOW 

VERY LOW 

LOW 

LOW 

LOW 

VERY HIGH 

HIGH 

HIGH 

HIGH 

VERY HIGH 

HIGH 

MEDIUM 

VERY_HIGH 

VERY HIGH 

HIGH 

LOW 

VERYJHIGH 

VERY HIGH 

MEDIUM 

HIGH 

HIGH 

VERY HIGH 

MEDIUM 

MEDIUM 

HIGH 

VERY HIGH 

MEDIUM 

LOW 

VERYJHIGH 

VERY HIGH 

Low 

HIGH 

LOW 

VERY HIGH 

LOW 

MEDIUM 

MEDIUM 

VERY HIGH 

LOW 

LOW 

MEDIUM 

HIGH 

HIGH 

HIGH 

HIGH 

HIGH 

HIGH 

MEDIUM 

HIGH 

HIGH 

HIGH 

LOW 

VERYJHIGH 

HIGH 

MEDIUM 

HIGH 

HIGH 

HIGH 

MEDIUM 

MEDIUM 

HIGH 

HIGH 

MEDIUM 

LOW 

VERYJHIGH 

HIGH 

LOW 

HIGH 

LOW 

HIGH 

LOW 

MEDIUM 

LOW 

HIGH 

LOW 

LOW 

MEDIUM 

TABLE4.10 RULES FOR SLAG CORROSION I 

RESISTANCE (2 of 2) 


CTI 

GHDT 

m 

W| 

MIXING 

TIME 

MIX 

TEMP 

PMT 

ST 

££ 

FP 

' t 1 ‘ f- , ' 

1 || 5 §§| 
lltjflll 

low 

high 

low 

high 

medium 

medium 

high 

medium 

low 

medium 

BSIKM 

low 

low 

low 

medium 

medium 

low 

very low 

high _ 

medium 

medium 

mm 

medium 

medium 

medium 

low 

hiqh 

low 

very low 

medium 

low 

low 

wSmmm 

low 

low 

high 

high 

medium 

medium 

medium 

medium 

medium 

high 


low 

medium 

high 

high 

high 

medium 

medium 

high 

low 

high 

WSM 

high 

high 

— . — . 

low 

J-Ua. 1 -! 

high 

’O'- — 

high 

medium 

medium 

medium 

medium 

high 


medium 

high 

low 

hiqh 

medium 

medium 

r 

medium 

high 

low 

medium 

WtESmm, 

high 

low 

low 

hiah 

high 

medium 

low 

low 

medium 

low 

medium 

medium 

low 

high 



high 

high 

high 

medium 

medium 

medium 

low 

■ r TT* _ 

very low 
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high 

medium 

medium 

j high 

low 

medium 

high 

medium 

low 

hinh 


. 

high 

medium 

J _ 

1 medium 

medium 

low 

low 

medium 

medium 

mjsaium 

high.— 

high 

: low 

1 _ low 

j medium 

high 

high 

low 

medium 

hiqh 

mprlh im 

high. , 

medium 

. J ow 

high 

medium 

medium 

medium 

medium 

medium 

4 — jj 1 1 

high 

mprftnm 

high_ 

__ m e d f u m 

Jow 

. high 

j medium 

high 

high 

medium 

medium 

hinh 


__hiflh . 

_ medium 

low 

‘ medium 

medium 

medium 

medium 

medium 

medium 

"‘U 1 1 

hinh 

meuium 

low 

low 

low 

] low __ 

1 . ^ ow 

j medium 

medium 

medium 

medium 

low 

hiqh 

hinh 

low 

low 

| medium 

j Sow 

medium 

medium 

low 

medium 

medium 

J — g 

high 

n| gl u 

low 

low _ 

' high 

j low 

medium 

low 

medium 

medium 

low 

hiqh 

medium 

low 

low 

medium 

J low 

medium 

low 

medium 

medium 

low 

hiqh 

low 

high 

medium 

low 

medium 

medium 

high 

high 

high 

low 

medium 

Vfirv Inuy 

high 

low 

low 

high 

low 

high 

high 

high 

medium 

medium 

' low 

high 

medium 

low 

low 

low 

medium 

high 

medium 

low 

medium 

low 

high 

medium 

low 

low 

low 

medium 

„ high 

medium 

low 

medium 

medium 

high 

low 

low 

medium 

medium 

low 

medium 

medium 

low 

low 

low 

medium 

low 

medium 

] medium 

medium 

medium 

medium 

low 

medium 

medium 

hiqh 

high 

low 

medium 

high 

low 

medium 

medium 

medium 

low 

medium 

; 

medium 

high 

low 

low 

high 

low 

low 

medium 

low 

low 

medium 

: v 

medium 

high 

low 

low 

medium 

low 

low 

medium 

medium 

low 

medium 

medium 

medium 

low 

high 

medium 

medium 

high 

medium 

high 

low 

high 

medium 

medium 

high 

low 

low 

high 

low 

veryjow 

low 

low 

medium 

medium 

high 

low 

low 

high 

high 

high 

medium 

high 

low 

medium 

veryjow 

high 

low 

low 

high 

medium 

high 

high 

high 

medium 

medium 

very low 

low 

low 

high j 

low 

medium 

medium 

high 

low 

low 

medium 

, '/ j 

low 

high 

-highj 

low 

low 

high 

high 

low 

medium 

medium 


low 

h!3h_. _ 

medium 

medium 

i low 

low 

medium 

medium 

low 

low 


low 

medium 

low 

low 

medium 

low 

low 

high 

low 

low 


medium 

...high 

low 

high 

medium 

medium 

medium 

low 

medium 

medium 


medium 

high 

low 

high 

medium 

medium 

medium 

low 

low 

medium 

high 

low 

medium 

low 

medium 

medium 

low 

low 

medium 

low 

medium 

1 

low 

low 

low 

medium 

medium 

low 

veryjow 

high 

low 

medium 

\ J 

low 

medium 

medium 

medium 

low 

high 

medium 

low 

medium 

medium ! 

\ ' 1 

low 

medium 

high 

high 

medium 

low 

medium 

high 

low 

medium J 

medium 

low 

medium 

jow 

high 

low 

medium 

medium 

high 

low 

medium t 

■ 

1 

low 

medium 

high 

low 

medium 

low 

low 

high 

low 

low 


low 

high 

low 

high 

medium 

medium 

medium 

medium 

medium 

medium | 

high 

medium 

medium 

high 

medium 

high 

medium 

high 

medium 

low 

low 

Veryjijgh 

medium 

low 

high 

medium 

medium 

medium 

medium 

medium 

low 

high | 


medium 

low 

medium 

medium 

medium 

low 

low 

low 

low 

low 

iSKSfiiS 

high 

medium 

medium 

medium 

high 

high 

high 

low 

medium 

low 


high 

high 

medium 

medium 

medium 

high 

high 

low 

medium 

low 


_ high 

high 

low 

medium 

medium 

medium 

medium 

medium 

medium 

medium 


. high 

high 

low 

high 

medium 

medium 

medium 

low 

medium 

medium 

, 

low 

low 

medium 

medium 

medium 

low 

low 

low 

low 

low 

m si 
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low 

low 

high 

J medium 

i medium 

i — - 

low 

medium 

medium 

low 

medium 

verv hinh 

low 

low 

— — 

1 medium 

1 medium 

4— - — 

medium 

low 

low 

high 

low 

high 

high 

high 

[ high 

!i 'gh 

. medium 

j medium 

— high 

high 

low 

medium 

hiqh 

verv hiah 

medium 

high 

high 

i . . * ow 

high 

low 

veryjow 

low 

low 

medium 

high 

low 

medium 

. h| gh 

low _ 

1 

j medium 

low 

low 

low 

low 

medium 

medium 

low 

medium 

high 

medium 

| medium 

low 

very low 

low 

low 

medium 

medium 

low 

medium 

i-hJSjL 

Sow 

medium 

low 

veryjow 

low 

low 

high 

hiqh 

low 

medium 

high 

Sow 

low 

medium 

medium 

medium 

low 

medium 

! — J 

low • 

low 

medium 

high 

| medium 

medium 

low 

medium 

medium 

low 

medium 

medium 

low 

high 

high.. . 

low 

medium 

low 

low 

high 

low 

medium 

hiqh 

high 

high _ 

high 

low 

medium 

medium 

medium 

medium 

medium 

high 

medium' 

high 

high 

high 

1 medium 

medium 

medium 

medium 

medium 

medium 

low 

hiqh 

high 

high 

medium 

| medium 

medium 

medium 

high 

low 

low 

high 

-J2-J J 

high 

high 

medium 

medium 

medium 

medium 

low 

low 

low 

low 

low 

very high 

medium 

low 

high 

low 

medium 

. .high , 

very_high 

low 

medium 

high 

medium 

medium 

low 

high 

medium 

medium 

low 

low 

low 

medium 

high 

very high 

medium | 

Jow 

high 

j medium j 

medium 

high 

low 

medium 

medium 

low 

medium 

medium 

low 

high 

j medium 

- high _ 

low 

medium 

low 

medium 

high 

, U. -,n. , 1. ill n. C„ u 1H, ,,i 

veryjow 

medium 

low : 

high 

L h'9 h - 

high 

high 

very_high 

low 

medium 

low 

medium 

low 

low 

high ; 

| high 

high 

high 

veryjtigh 

low 

medium 

medium 

medium 

low 

L !0W | 

high | 

high 

high 

high 

very_high 

medium 

low 

medium 

low 

medium 
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low 

medium 

- high 
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low 

low 

medium 

high 

medium 

low 

medium 

medium 

medium 

low 

low 

medium 

medium 


high 

medium 

i°w 

medium 

low 

medium 

medium 

medium 

medium 

medium 

medium 

-high 

medium 
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medium 

low 

medium 

medium 

low 

low 

high 

low 

medium 

medium ^ 

low 

.. high... 

low 

medium 

medium 

low 

low 

medium 
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Table 4.1 1 Actual Data of 86 sets Obtained from RSP Converted to Very High, 
High, Medium, low , very low (3of 3) 


APPENDIXB 


RULE BASE FOR EXPERT BASED DATA OttTAmirn mN G Tm 

PROGRAM ~ 


ghdt = high and rt -- high and 

cti = high ==> class = very high 

cti = medium ==> class = very high 
cti = low ==> class = high 

il = medium and cti = high ==> class = high 

cti = medium ==> class = high 
cti = low ==> class = medium 

rt " low and 

cti = high ==> class = medium 
cti - medium ==> class = medium 
cti = low ==> class = low 

ghdt medium and rt ~ high and 

cti = high ==> class = high 
cti = medium ==> class = high 
cti = low ==> class = medium 

rt medium = 

~> class = medium 

rt = low and 

cti = high ==> class = medium 
cti = medium ==> class = low 
cti = low ==> class = low 

ghdt = low and rt = high and 

cti - high ==> class = medium 
cti = medium ==> class = medium 
cti = low — > class = low 

rt medium and 

cti - high ==> class = low 
cti = medium ==> class = low 
cti = low — > class = veryjow 

i! 

o 

II 

II 

V 

class = very_low 


Figure 4.5: Decision Tree for Grain Temperature of Table 4.2 
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gt = very high > class very high 

gt = high and mixer temp high “> class = high 

mixer temp = medium class = high 

mixer temp - low and mixjime - high => class = high 

mixjime = medium ==> class = high 
mixjime = low ==> class = medium 

gt = medium and mixer temp high — > class = medium 

mixer temp medium ™> class = medium 

mixer temp low and mixjime = high — > class = low 

mixjime = medium => class = low 
mixjime = low ==> class = medium 


gt - low and mixer temp high and mixjime - high — > class = medium 

mix time - medium ==> class = low 
mix time = low ~=> class = low 

mixer temp medium class = low 
mixer temp low : > class - low 

gt = veryjow and mixer Jemp - high and mixjime - high ==> class = low 

mixjime = medium ==> class = low 
mixjime = low — > class = veryjow 

mixer temp - medium and mixjime = high => class = low 

mixjime = medium ==> class = 

mix time = low ==> class = 


mixer temp low —> class * veryjow 



low 

very low 
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mix temp very high and rolls high and st = high ==> class = low 

st = medium ==> class = low 
st - low ==> class = medium 

rolls - medium and st = high ==> class = medium 

st = medium — > class = medium 
st = low ==> class = high 

rolls = low and st = high — > class = high 

st = medium ==> class = high 
st = low ==> class = very_high 

mix temp high and rolls = high ==> class = medium 

rolls medium and st - high ===> class = medium 

st - medium ==> class = medium 
st “ low ==> class = high 

rolls low -^> class - high 

mix temp medium and rolls high r ~> class ~ low 

rolls medium “> class ® low 
rolls low “> class = medium 

mixjemp = low class = low 

mixjemp = veryjow and rolls = high ==> class = veryjow 

rolls -medium ==> class = veryjow 
rolls - low and st = high ==> class = veryjow 
st = medium — > class = low 
st = low ==> class = low 


Figure 4.7: Decision Tree for Press Mix Tempera ture of table 4.4 




graphite ver> low and pint high and fp = high and pitch = high => class = veryjhigh 

pitch = medium ==> class = veryjhigh 
pitch = low => class = high 

fp - medium and pitch = high ==> class = veryjhigh 

pitch = medium — > class = very_high 
pitch = low — > class = high 

fp = low and pitch = high ==> class = high 

pitch = medium ==> class = high 
pitch = low ==> class = medium 

pmt =■ medium and fp = high — > class = high 

fp = medium and pitch = high ==> class = high 

pitch = medium ==> class = high 
pitch = low ==> class = medium 

fp “ low and pitch = high ===> class = high 

pitch = medium ==> class = medium 
pitch = low ==> class = medium 

pmt low and pitch - high “> class = medium 

pitch medium and fp = high — > class = high 

fp = medium — > class = medium 
fp = low ==> class = medium 

graphite very high and fp high and pmt : = high and pitch = high ==> class = low 

pitch = medium ==> class = low 
pitch = low ==> class = very Jow 

pmt = medium and pitch = high ==> class = low 

pitch = medium ==> class = veryjow 
pitch = low ==> class = veryjow 

pint « low --> class = very_low 

fp medium and pmt = high and pitch = high ==> class = low 

pitch = medium ==> class = veryjow 
pitch = low — > class = very_low 

pmt = medium ==> class = veryjow 

pmt = low ==> class = veryjow 

fp - low ™> class = very__Iow 



graphite high and 

ip high and 

pmt - high and pitch = high ==> class = medium 
pitch — medium ==> class = low 
pitch = low ==> class = low 



pmt = medium ==> class = low 
pmt = low ==> class = low 


fp = medium and pmt = high ==> class = low 

pmt = medium ==> class = low 
pmt = low ==> class = very low 


fp = low and 

pitch = high and pmt = high ==> class = low 

pmt = medium ==> class = low 
pmt = low ==> class = very low 



pitch = medium ==> class = veryjow 
pitch = low ==> class = very_Iow 

graphite medium and pmt high and fp 

tp= 

high ==> class = medium 
medium and pitch = high ==> class = medium 

pitch = medium ==> class = medium 
pitch = low ==> class = low 


ip = 

low and pitch = high ==> class = medium 
pitch = medium ==> class = low 
pitch = low ==> class = low 


pmt = medium and 

fp = high ==> c i ass = medium 

fp = medium and pitch = high ==> class = medium 

pitch = medium ==> class = medium 
pitch = low ==> class = low 



fp = low and pitch = high ==> class = medium 

pitch = medium ==> class = low 
pitch = low ==> class = low 


pmt = low and 

fp = high and pitch = high => class = low 

pitch = medium ==> class = medium 
pitch = low ==> class = low 



fp = medium ==> class = low 

fp = low and pitch = high — > class = low 

pitch = medium ===> class = low 
pitch = low ==> class = veryjow 

Figure 4.8: 

Decision Tree for Green Bulk Density of table 4.5(2 of 3) 
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graphite low and pitch high and 

pint = high and fp = high ==> class = high 

fp = medium ==> class = high 
fp = low ==> class = medium 


pint = medium ==> class = high 

pint = low and fp = high => class = medium 

fp = medium ==> class = medium 
fp = low ==> class = low 

pitch = medium and 

pmt = high and fp = high ==> class = high 

fp = medium ==> class = high 
fp = low ==> class = medium 


pmt - medium — > class = medium 
pmt = low ==> class = medium 

pitch low and 

pmt = high and fp = high ==> class = medium 

fp = medium ==> class = medium 
fp = low ==> class = low 


pmt “ medium ==> class = medium 

pmt = low and fp = high ==> class = medium 

fp = medium ==> class = medium 
fp = low ==> class = low 


Figure 4.8: Decision Tree for Green Bulk Density of table 4.5( 3 of 3) 





gbd rr medium and pitch high and tempering = perfect ==> class = high 

tempering = under => class = high 
tempering = over ==> class = medium 

pitch medium and tempering = perfect — > class = medium 
tempering = under ==> class = medium 
tempering = over ==> class = low 

pitch low and tempering = perfect ==> class = low 

tempering = under ==> class = medium 
tempering = over ==> class = low 

gbd = low and pitch - high and tempering = perfect ==> class = medium 

tempering - under ==> class = medium 
tempering = over — > class = low 

pitch = medium and tempering = perfect ==> class = medium 
tempering = under ==> class = medium 
tempering = over ==> class = low 

pitch low and tempering = perfect ==> class = low 

tempering = under ==> class = low 
tempering = over ==> class = veryjow 

gbd = very low and tempering perfect ~> class = veryjow 

tempering under and pitch = high — > class = low 

pitch = medium ==> class = low 
pitch = low — > class = veryjow 

tempering = over =-> class = very_low 

gbd = very high and pitch = high and tempering = perfect ==> class = very_high 

tempering = under ==> class = veryjiigh 
tempering = over ==> class = high 

pitch ~ medium and tempering = perfect => class = veryjhigh 
tempering = under ==> class = high 
tempering = over ==> class = high 

pitch low and tempering = perfect =-> class = high 
tempering = under ==> class = medium 
tempering = over ==> class = very_high 

gbd = high and pitch - high and tempering - perfect => class = very_high 

tempering = under ==> class = high 
tempering = over — > class = high 

pitch = medium and tempering = perfect — > class = high 
tempering = under ==> class = medium 
tempering ~ over class = medium 

pitch - low “> class = medium 


Figure 4.9: Decision Tree for Cold Cr ushing Strength of table 46 
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gbd ; 

v er ) ^igh an< ^ pitch - high — > class = medium 

pitch - medium — > class = low 
pitch = low ==> class = tow 

ghd ^ 

high and pitch = high ==> class = high 

pitch = medium ==> class = medium 
pitch = low ==> class = low 

gbd 

medium and pitch = high ==> class = high 

pitch = medium ==> class = medium 
pitch = low ™> class = medium 

gbd 7 

low and pitch = high — > class = high 

pitch = medium ==> class = medium 
pitch “ low — > class = medium 

gbd 

very low -■"> class = high 


Figure 4.10: Decision Tree for Coked Porosity of table 4.7 





high and ccs high and gbd very low and 


ic very__high — > class = vvjiigh 
rc = high ==> class = veryjhigh 
rc = medium ==> class = high 


gbd high --> class = very_high 

gbd medium and rc = high ==> class = vvjiigh 


gbd ^ low and 


rc - medium — > class = vvjiigh 
rc = low ==> class = very_high 

rc = veryjhigh ==> class = vvjiigh 
rc = high ==> class = vvjiigh 
rc = medium ==> class = very_high 
rc = low — > class = very__high 


gbd = veryjiigh ==> class = high 


ccs medium and gbd very low and 


rc = veryjiigh ==> class = veryjiigh 
rc = high ==> class = veryjiigh 
rc = medium ==> class = medium 


gbd high ==> class = veryjiigh 

gbd - medium and rc = high => class = vvjiigh 

rc = medium ==> class = vvjhigh 
rc = low ==> class = veryjiigh 

gbd low and rc - veryjhigh ==> class = vv_high 
rc - high ==> class = vvjiigh 
rc = medium ==> class = veryjiigh 
rc = low ==> class = very_high 

gbd - veryjhigh ==> class = high 

ccs = low and gbd - veryjow and rc = veryjiigh ==> class = veryjiigh 

rc = high ==> class = very_high 
rc ~ medium — > class = medium 

gbd =* high =>-> class = high 

gbd - medium and rc = high ==> class = vvjiigh 

rc = medium ==> class = vvjiigh 
rc = low ==> class = very_high 

gbd = low and rc = veryjiigh — > class = vv__high 
rc = high ==> class = veryjiigh 
rc = medium ==> class = veryjhigh 
rc = low ==> class = high 

gbd = very__high ==> class - medium 


Figure 4.11: Decision Tree for HMOR of table 4.8 (1 of 3 
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medium and ecs high and gbd - very low and rc = veryjiigh ==> class = high 

rc = high ==> class = high 
rc = medium ==> class = medium 

gbd high "> class = high 

gbd medium and rc = high ==> class - vvjiigh 

rc = medium => class = vvjiigh 
rc = low ==> class = high 

gbd = low and rc = veryjiigh ==> class = vv_high 
rc = high ==> class = veryjiigh 
rc = medium ==> class = very_high 
rc = low ==> class = high 


gbd = veryjhigh — > class = medium 

ccs medium and gbd = very low and rc = veryjiigh ==> class = high 

rc = high ==> class = high 
rc = medium ==> class = medium 

gbd = high class = high 

gbd = medium and rc = high ==> class = vvjiigh 

rc = medium ==> class = veryjhigh 
rc = low ==> class = high 


gbd low and rc = veryjiigh ==> class = vv__high 
rc = high ==> class = very_high 
rc - medium ==> class = veryjiigh 
rc “ low ==> class = high 

gbd “ veryjhigh ™> class = medium 

ccs = low and gbd = veryjow and rc = veryjiigh ==> class = high 

rc = high ==> class = medium 
rc = medium ==> class = low 


gbd ~ high ==> class = high 

gbd - medium and rc = high ==> class = very_high 
rc = medium ==> class = high 
rc = low ==> class = high 

gbd =* low and rc = veryjiigh ==> class = veryjhigh 
rc = high => class = high 
rc = medium ==> class = high 
rc = low — > class = medium 


gbd - veryjiigh — > class = medium 


Figure 4.11: Decision Tree for HMOR of ta ble 4.8 (2 of 3) 
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low and rc \ er\ high and gbd low and 


ccs ~ high ==> class = high 
ccs = medium => class = high 
ccs = low — > class = medium 


gbd r: veryjow class = medium 

rc high and ccs = high and gbd - medium ==> class = high 

gbd = low ==> class = high 

gbd = veryjow — > class = medium 


ccs - medium ==> class = medium 
ccs = low -=> class = medium 

rc medium and gbd ~ medium and ccs = high ==> class = medium 

ccs = medium ==> class = medium 
ccs = low ==> class = low 

gbd low and ccs « high ==> class = medium 

ccs = medium ==> class = medium 
ccs ~ low ==> class = low 

gbd very low and ccs = high — > class = low 

ccs = medium ==> class = veryjow 
ccs = low ==> class = veryjow 


rc = low and gbd ~ high and ccs = high ~> class = medium 

ccs = medium ==> class = medium 
ccs = low ==> class = low 

gbd - medium and ccs = high ==> class = medium 
ccs = medium => class = low 
ccs = low — > class == low 

gbd - low and ccs = high ==> class = low 

ccs medium ==> class = veryjow 
ccs « low ==> class = veryjow 


rc = very low == 5 > class = low 


Figure 4,11: Decision Tree for HMOR of table 4.8(3 of 3) 
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chp ~ high and cp 

hipli and re medium ==> class = medium 
re low => class - low 
re - very Jaw => class = low 
rc“ veryjiigh — > class = high 
re high ==> class = high 

cp 

medium and re medium =~> class ~ high 
re low class = medium 

1C ^ er\ low ==> class - medium 
re - verv high ~> class « veryjiigh 
re = high ~=> class = high 

cp 

- low and re - medium => class ~ high 
re r low ~=> class = medium 
rc = veryjow class = medium 

re = ven high =-> class - very high 
re t: high -=> class = veryjiigh 

chp = medium and 

r ■ medium and cp ~ high ==> class = medium 

cp - medium —> class = medium 
cp - low =-> class - high 

re 

low and cp - high ™> class - low 

cp medium -=> class = medium 
cp low ~=> class = medium 

IX 

ven low and cp high -=> class = veryjow 
cp medium “-> class = low 
cp ~ low — > class = low 

re 

ven high and cp high —> class = high 

cp - medium =~> class = high 
cp - low class = veryjiigh 

re 

high and cp = high => class = high 

cp “ medium —> class = high 
cp - low — > class = veryjiigh 

chp = low and re - 
re = 

medium -=> class = low 

low and cp - high --> class = veryjow 

cp - medium — > class - veryjow 
cp low class = low 

re 

ven low and cp high -=> class = veryjow 

cp = medium “=> class - veryjow 
cp " low =*-> class - low 

re 

very high and ep ~ high ~> class - low 

cp = medium => class = medium 
cp = low => class = medium 

re = 

high and cp = high — > class ~ low 

cp = medium => class = low 
cp = low — > class = medium 
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rc ven low and mp high and gbd = very Jiigh -=> class = medium 

gbd - high ==> class = low 
gbd = medium => class = low 
gbd = low — > class = low 

mp medium and gbd = veryjiigh => class = low 
gbd - high ==> class = veryjow 
gbd - medium ==> class = low 
gbd = low =» class = veryjow 

mp low > class = veryjow 

rc = low and gbd high and mp = high =» class = medium 

mp = medium ==> class = low 
mp ~ low ==> class = low 

gbd medium and mp = high ==> class = medium 
mp ~ medium ==> class = low 
mp = low ==> class = low 

gbd low and mp = high =-> class = low 

mp - medium ==> class = veryjow 
mp - low ==> class = veryjow 

rc high and gbd medium and mp = high ==> class = veryjiigh 

mp - medium ==> class = high 
mp - low class = medium 

gbd low and mp ~ high —> class - high 

mp - medium =-> class = high 
mp “ low => class = medium 

gbd very low and mp ~ high ==> class = medium 

mp = medium ==> class = medium 
mp = low ==> class = low 

rc : = medium and gbd - medium and mp = high => class = high 

mp = medium ==> class = high 
mp = low ==> class = medium 

gbd = low and mp ~ high ==> class = medium 
mp = medium ==> class = low 
mp * low ==> class = low 

gbd veryjow and mp = high ==> class « medium 
mp =* medium ==> class = low 
mp - low ==> class ® low 

rc ~ very high and gbd • low and mp - high ~~> class = veryjiigh 

mp® medium —> class » veryjiigh 
mp = low ==> class = high 

gbd = veryjow and mp = high ==> class ® very_high 
mp = medium === > class = high 
mp - low ==> class = medium 
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Fig 4.14 i rcc for data s et of 86 samples shown on two pagesd of 2^ 


MTM = H & PN1 1 I ( I \SS M 

PMT-'H & GHDT M & Ml M & FP = 1. & CTI = M ==> CLASS = VH 

CT1 = H —> CLASS = M 

PMT - VI & CTI M ii GHDT -M -■ > CLASS = VL 

0111,1 11 & MT = L & MIX “ L & ST = L & ROLLS = L & FP = M & RT = L ==>CLASS = M 

PMT M A (.111)1 M -( I ASS - M RT = H => CLASS = H 

( .111)1 II ( I ASS M 

i .111) I 1 CLASS -- VL 

PMT = Vi I A si | ( LASS M 

S I M l I ASS | 

MTM = M & MIX = 1 & ROLLS I & CTI - M -=> CLASS * M 

( 'l l L & RT = L & GHDT = M & MT = L => CLASS = VL 

MT = M => CLASS = H 
GHDT = L ==> CLASS = M 

RT = M & MT - L ==> CLASS = L 
MT = M — > CLASS = H 

RT li & GHDT = L & PMT = M & ST = M & MT = L => CLASS = M 

MT = M => CLASS = VH 
GHDT = M & FP = L ==> CLASS = M 
FP = H => CLASS = H 
FP = M => CLASS = M 
GHDT = H ==> CLASS = H 

("II 1 1 & MT M & FP = I. & GHDT = L ==> CLASS = L 

GHDT = M ==> CLASS = VH 

ROl.l S M A GHDT I. A MT = M & CTI = L -o CLASS -L 

CTI = M ==> CLASS = VH 

MIX = 11 & FP = M » CLASS = VI. 

FP = II & ROLLS - I. ==> CLASS = M 

ROLLS = M & CTI = IT & RT = L ==> CLASS = M 
RT = H => CLASS = VH 
CTI = M ==> CLASS = M 

FP I . & CTI = M &. MT = M & ROLLS = M & GHDT = L => CLASS = M 

GHDT = M ==> CLASS = H 

("IT II * CLASS - M 

MIX- M & GHDT II & ROLLS M & PMT * L —» CLASS **M 

PMT = M & RT = L & MT = M ==> CLASS = H 

MT = H & ST = M ==> CLASS = H 

ST = L & FP = M & CTI = M => CLASS = M 
CTI = H => CLASS = H 

RT = H & CTI = H & ST = M & MT= L =>CLASS = M 

MT = M => CLASS = H 

ROLLS = L & ST = M => CLASS = M 
ST = H => CLASS = L 
ST = L => CLASS = H 

GHDT = L & ST = I. ==> CLASS = H 

ST = M & PMT = L ==> CLASS = L 
PMT = M ==> CLASS = H 


GHDT = M & ROLLS = L ==> CLASS = H 

ROLLS = M & CTI = M & RT = M => CLASS -M 
RT = L => CLASS = VL 
CTI = H & MT = H => CLASS = M 


MT = M & RT = L & PMT = M => CLASS -L 


n\ at ~ r 


n T A CQ = N/f 




UK 4r , H ll? c f° r d . f*t a set of 86 sam p les shown on two pages(2 of 2) 


MTM = l & KOI 1 S M ft I'M I II & MIX 11 & ] )• M & CTI = H => CLASS = L 

CTI L ==> CLASS = H 

I ’M i M Cl ASS - M 

ROLLS I ft I’M L II & CTI H & RT = M ==>CLASS = M 

RT = L & GHDT = M & MT = L & MIX = M & ST = M & FP = M =>CLASS = M 

FP = M => CLASS = L 

CTI M ==> CLASS = H 

PM I M & MIX - 1. — > CLASS = M 

MIX M & GHDT = L ==> CLASS = M 

GHDT = M & RT = H => CLASS = L 

RT = L & MT = M => CLASS = L 
MT = IT => CLASS = H 
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Appendix C 


PLANT DATA FOR PROCESS UP TO GBP 


CTI 

GHDT 

RT 

MIXER i MIXING 
TEMP 1 TIME 

MIX 

TEMP 

PMT 

ST 

NO OF 
ROLL 

FP 

GBD 

414 

519 

6 

126 16 

153 

132 

85 

1 

1484 

3.09 

385 

532 

5 

128 15 

135 

128 

90 

1 

1353 

3.07 

380 

480 

6 

123 

16 

140 

130 

28 

1 

1345 

3.11 

361 

494 

6 

126 

16 

135 

128 

108 

0 

1345 

3.05 

403 

511 

6 

120 

14 

145 

138 

45 

0 

1484 

3.06 

364 

502 

6 

124 

14 

125 

121 

20 

0 

1345 

3.05 

353 

454 

6 

120 

14 

125 

120 

20 

0 

1484 

3.07 

365 

466 

6 

1 16 

14 

140 

125 

60 

0 

1396 

3.05 

360 

41 1 

7 

123 

14 

144 

144 

10 

1 

1476 

3.09 

365 

381 

9 

120 

12 

137 

130 

50 

1 

1360 

3.06 

368 

369 

7 

1 18 

14 

142 

130 

85 

0 

1484 

3.06 

380 

490 

7 

105 

16 

135 

135 

20 

1 

1484 

3.08 

364 

492 

6 

135 

15 

153 

140 

40 

0 

1491 

3.07 

359 

484 

7 

133 

16 

143 

135 

210 

0 

1462 

3.12 

404 

505 

6 

127 

16 

136 

135 

12 

0 

1258 

3.14 

374 

521 

6 

129 

17 

135 

134 

19 

0 

1353 

3.1 

354 

500 

6 

129 

17 

132 

131 

12 

0 

1360 

3.13 

364 

539 

6 

129 

16 

134 

132 

28 

0 

1338 

3.12 

359 

539 

6 

131 

14 

138 

135 

31 

0 

1353 

3.12 

373 

587 

6 

129 

14 

132 

128 

128 

1 

1345 

3.12 

460 

380 

7 

125 

16 

145 

135 

30 

1 

1491 

3.12 

357 

466 

8 

131 

16 

146 

136 

38 

1 

1484 

3.06 

349 

475 

8 

137 

16 

146 

138 

83 

1 

1476 

3.11 

346 

472 

8 

135 

20 

145 

140 

90 

2 

1484 

3.05 

382 

528 

6 

118 

18 

160 

140 

45 

1 

1476 

3.15 

376 

498 

4 

125 

17 

129 

125 

5 

1 

1345 

3.08 

395 

490 

7 

123 

28 

135 

130 

5 

1 

1353 

3.08 

377 

510 

7 

129.6 

26 

140 

131 

20 

1 

1345 

3.04 

363 

479 

7 

135 

27 

142 

132 

30 

1 

1345 

3.08 

_ 355 

477 

7 

136 

27 

140 

130 

40 

1 

1353 

3.06 

. 330 

614 

5 

120 

7 

120 

115 

76 

1 

1411 

3.1 


Table 4.12 31 Data Samnle Points of RSP Data 
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GHD'I 
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'"551.5 
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_420.I 
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435.8 
428 
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“To 1 
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MIXER 

TEMP 
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MIXING 

TIME 

8’ 

8 
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~ 11 ~ 


M 

7 I 


5 

'4 
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MIX 

TEMP 
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130 
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PMT ST 
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112 30 

122 35 

120 175 

122 48 

120 108 
115 20 

122 65 

128 31 

115 26 

128 20 
125 30 

127 45 

123 45 

' 122 33 

117 42 

120 43 

120 38 

128 115 
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125 35 

123 25 
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125 40 
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112 20 
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116 80 
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118 40 
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NO OF 
ROLLS FP 
0 195 


FP 
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1527.27 
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1410.91 
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1541.82 
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1410.91 
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1403.64 
1403.64 
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1410.91 
1432.73 
1410.91 
1425.45 
1418.18 
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0 
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441 

8 

119 

18 

135 

130 

37 

0 

185 

1345.45 

3.17 

321.2 

425.1 

8 

124 

8 

130 

125 

31 

0 

218 

1585.45 

3.15 

'330 

408 

6 

122 

12 

116 

114 

20 

0 

188 

1367.27 

3.12 

Ho - 

520 

6 

118 

16 

140 

128 

13 

1 

192 

1396.36 

3.11 

417 

602 

6 

118 

9 

143 

130 

25 

1 

192 

1396.36 

3.1 

401 

580 

5 

124 

8 

134 

120 
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oo 

VO 
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1 

196 

1425.45 

3.14 
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8 
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0 
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3.16 
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11 
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0 
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10 
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78 

0 

202 

1469.09 

3.15 
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8 

124 

12 
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25 

1 
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14 
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0 
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3.15 
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482 
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10 
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6 

0 
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9 
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13 
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12 
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29 

0 
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3.08 
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9 
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8 
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1 
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8 
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1 
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0 

204 
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0 
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418 

8 
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204 
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11 
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17 
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120 
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323 

401 
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1 
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1389.09 

3.09 
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392 
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148 

142 
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1 

199 
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0 
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392 
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1 
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3.14 
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467 
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8 
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85 

1 
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1381.82 

3.09 

376 

505 

5 

117 

8 

135 

124 

72 

1 

190 

1381.82 

3.03 

356 

487 

6 

123 

8 

138 

125 

79 

1 

190 

1381.82 

3.13 

368 

494 

5 

125 

10 

136 

124 

72 

0 

190 

1381.82 

3.14 

377 

507 

6 

127 

7 

140 

128 

91 

0 

190 

1381.82 

3.13 

390 

501 

5 

125 

9 

142 

129 

62 

0 

190 

1381.82 

3.12 

383 

583 

5 

128 

8 

138 

124 

25 

0 

190 

1381.82 

3.12 

391 

506 

5 

122 

11 

130 

118 

15 

1 

195 

1418.18 

3.1 

394 

522 

4 

124 

5 

137 

121 

30 

1 

198 

1440 

3.11 

389 

516 

5 

121 

6 

127 

125 

5 

0 

203 

1476.36 

3.06 

_ 375 

501 

4 

126 

4 

128 

126 

5 

0 

196 

1425.45 

3.13 


Table T4.13 1 Actual Data of 86 sets Obtained from RSP 



Appendix D 


LISTING OF ALL RULES AUTOMATICALLY GENERATED BY CART 


treekrilks 

1. ins I • 5X.5 ( il IDT <- 537.5 && CTI <= 342 && MIXING TIME <= 10.5 && 

MIX I R I EMPERATURE <= 1 14.6) {Class = low} 

2. If (ST • = 58.5 &&GHDT <= 537.5 &&CTI <= 342 &&MIXING TIME <= 10.5 && 

MIXER TEMPERATURE > 1 14.6) {Class = high} 

3. If (ST <= 58.5 &&GHDT <= 537.5 &&CTI <= 342 &&MIXING TIME > 10.5) 

I Class very high] 

4. If (ST • 58.5 && CTI > 342 &&GHDT <= 532 && FP <= 1337.5){Class = very high} 

5. If (ST - 58.5 &&C1 1 > 342 &&GHDT <= 532 &&FP > 1337.5 &&PMT <= 121.5) 

{Class very low} 

6. If(GIII>l • 532 &&FP > 1337.5&&PMT > 121.5 &&CTI > 342&&CTI <= 405.5&& 

ST • 29 &&M1XER TEMPERATURE <= 129.3 &&MIX_TEMP <= 137.5) 

{Class low} 

7. If ((SHOT <= 532 &&FP >1337.5 &&PMT> 121.5 &&CTI > 342 &&CTI <=405.5 && 

ST <= 29 &&MIXER TEMPERATURE <= 129.3 &&MIX_TEMP > 137.5) 

{Class = medium} 

8. If (GIIDT < 532 &&FP > 1337.5 &&PMT> 121.5 &&CTI > 342 &&CTI<=405.5 && 

ST <= 29 &&MIXER TEMPERATURE > 129.3) {Class = very low} 

9. If (FP > 1337.5 &&PMT > 121.5 &&CTI > 342 &&CTI<= 405.5 &&ST > 29 &&ST 
< 58.5 && CrI IDT <- 522 &&MIXING TIME <= 6) {Class = medium} 

10. If (FP > 1337.5 &&PMT > 121.5 &&CTI > 342 &&CTI <= 405.5 &&ST > 29 &&ST 
<= 58.5 &&GI IDT <= 522 &&MIXING TIME > 6) {Class = low} 

11. If (FP > 1337.5 &&PMT > 121.5 &&CTI > 342 &&CTI <= 405.5 &&ST>29&&ST 
<= 58.5 &&GHDT > 522 &&GHDT <= 532) {Class = high} 

12. If (ST <= 58.5 &&GHDT <= 532 &&FP > 1337.5 &&PMT> 121.5 &&CTI > 405.5) 

{Class = medium} 

13. If (ST <= 58.5 &&CTI > 342 &&GHDT > 532 &&GHDT <= 537.5){Class - very high} 
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14. IflS'l • 58,5 &&Ci I • 409 &&PMT <= 130.5 &&GHDT>537.5 &&GHDT <=81. 5) 

\ Class highj 

15 _ If { S I • 58. 5 &&CTI <- 409 &&PMT <= 130.5 &&GHDT > 581. 5){Class = medium} 

16. IffS l • 58.5 &&C il 10 1 > 537.5 &&CTI <= 409 &&PMT > 130.5){Class = medium} 

17. IftS l ' 58.5 &&GI IDT > 537.5 && CTI > 409) {Class = medium} 

18. If (SI 58.5 &&GHDT <= 413.55) {Class = low} 

19. If (ST • 58.5 &&GHDT > 413.55 &&GHDT <= 505.8 &&CTI <= 289.5) 

| Class ; medium*, 

20. If (ST - 58.5 &&GHDT> 41 3.55 &&GHDT<= 505.8 &&CTI> 289.5) 

; Class very low} 

21. If(S I • 58.5 && (II IDT > 505.8) {Class = high} 



Appendix E 


Ol’TPrTS OBTAINED USING CART ANALYSTS 


Variable Importance 

The variable importance as analyzed by Cart is shown below, 


VARIABLE % 

IMPORTANCE 

r 

( T 1 1 1 ) 1 

100,000 

r 

MIX IEMP 

89.650 

r 

("11 

77.924 

r 

MIXER 1 IMP 

60.601 

r 

ST 

55.449 

r 

PM I 

54.715 

r 

MIXING 1 1MI 

54.137 

r 

RT 

47.867 


IP 

30.409 

'r 

NO OF ROLE 

11.518 


However as per the expert opinion the order of variable importance should have 

been as under, 

^ PMT 

> MIX TEMP 

> MIXING TIME 

> CTI 

^ GIIDT 
'r RT 
^ ST 

> EP 

> NO of ROLLS 
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1 uc ydm L IUU ,c P Icscm ine P^centage of a class in the population as it increases 
from 0 % to 100 %. I his is useful in the sense that one gets the overall picture of how 
many instances ol each class gets included as the population increases to 100%. Here 

Class 1 to class 4 represents GBD Very Low, Low, Medium, High 


Gain* Chart 


T erminal N 


Chart 


T erminal N 


Terminal N 


T erminal Ni 


% Population 


% Population 








1 or a total ni . I terminal Nodes generated by the Cart tree we have 21 Rules sets 
one ior each terminal node I he % ( lass composition of these nodes based on these 21 
rules is dispkned by Bar Charts and Pie charts below for Node 1 to Node 6. Here the 

classes ! to 6 ( VI . 1 . M. 1 1. V! 1 1 are represented by colour codes. 







Node 3 





Appendix F 


PROBABILTY PREDICTION FOR GBP TAKING A SPECIFIC 


EXAMPLE 





ijvitflipti 1 


Console 


Cftafi 


^ G ram_h e ate r_d ru m_te m p 

R?fntlon_tim« cSusion.chamberJemp 

__ 

J§r«mp 


No" of rolls Mixjttmp Storage Jimt 


f ; Prass_mixj«mp { .) 

Formin0j>ressur«r Pitch 


jJLirz: . 

fttttvariafel* i Sait Function | Edit Network 


Posterior distribution: 

probability ( "Grain Jemp" ) { //I variable(s) and 5 values 
table 

0.0 // p(Very_high | evidence ) 

0.0 //p(High| evidence) 

1.0 //p(Medium | evidence) 

0.0 //p(Low | evidence) 

0.0; //pC/ery low | evidence); 

} 

To observe a node, click on it. 

To query on a particular node, click on it. 


Posterior distribution: 
probability ( M Mixer_temp" ) { , 
table 

1.0 

0.0 

o.o; 

} 

Posterior distribution: 
probability ("Mixing Jme"){ 
table 

0.0 

1.0 


//I variable(s) and 3 values 

// p(High | evidence ) 
//p(Medium| evidence) 
//p(Low | evidence); 


//I variable(s) and 3 values 

//p(Highl evidence) 
//p(Medium | evidence) 
//p(Low | evidence); 


Posterior distribution: 
probability ( "Mix Jemp" ) { //I 
table 


variable(s) and 5 values 

// p(Very_high I evidence ) 
//p(High | evidence) 
//p(Medium | evidence) 
//p(Low | evidence) 
//p(VeryJow | evidence); 


Ftf» 6.6 Probability Prediction Obtained for Mix Temp (Medium 1) 
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m 

Gram_h#«*r^drum^t#mp 
R “n«lon_thn* cSu„on_ch3mb«r_t« mp 


MIxinsLtim 


3f 

m Grimj 




No^of^roHs Mix_t*mp 5tor*g*_tim* 

m j Prii^mix^Ump {..■>} 

Forming jpr*%sur+ Pitch 

3 k „ jl 


Posterior distribution ““ ’ “ “ — ,m fflf 11 r ^ 

probability ( "Mixjemp" ) { m variable(s) and 5 values 
table 

0,0 H p(VeryJiigh [ evidence ) 

0 0 tfp(High | evidence ) 

1.0 tf p(Medium | evidence ) 

0.0 //p(Low | evidence) 

j 0.0; //p(Veryjow | evidence); 

To observe a node, click on it. 

To query on a particular node, click on it. 

Posterior distribution: 

' probability ( "Storagejime" ) { //I variable(s) and 3 values 
table 

0.0 //p(High | evidence ) 

0.0 //p(Medium | evidence ) 

! 1-0; //p(Low | evidence); 

Posterior distribution: 

probability ( "No„or_rolls'* ) { m variable(s) and 3 values 
table 

0.0 //p(High | evidence) 

0.0 //p(Medium | evidence) 

1.0; //p(Low | evidence); 

) 

Posterior distribution: 

probability ( "Pre$s_mix_temp" ) { //I variable(s) and 3 values 
table 

0.0 // p(High | evidence ) 

1.0 //p(Medium | evidence) 

0.0; //p(Low | evidence); 


««»Witf*6»# 1 ftMtryndNMt | ftdrtftftfaHft 


Fig 6.7 Probability Prediction Obtained for PMT (Medium=l) 


m 

Grim^h**t#r - drum t#mp 


PI |» 1 1"^ I Pu> .laviBivtsCo^n 


Filff Options Help 

j To query on a particular node, click on it 


Retention time 


Combusionj:hamber_temp 


Mlxmg_time GriTn_t#mp 

• e 

No_of rolls Wix^temp 


Mix«rjt«mp 


Storage^tlme 


Formlng a j>r**fur# 


Fresi^mlx^temp 


Green bd 


Graphite 


Posterior distribution: 

probability ( "Pitch" ){ //I variable(s) and 3 values 
table 

| 0.0 //p(High | evidence) 

1.0 ffp(Medium I evidence) 

0.0; //p(Low | evidence); 

> 

Posterior distribution: 

probability ( "Graphite" ) { //I variable(s) and 5 values 
table 

0.0 // p(Very_high | evidence ) 

0.0 //p(High | evidence) 

1.0 // p(Medium | evidence ) 

0.0 //p(Low | evidence) 

0.0; // p(Veiy_low | evidence ); 

} 

Posterior distribution: 

probability ( "Formingjaressure" ) { //I variable(s) and 3 values 
table 

0.0 // p(High | evidence ) 

1.0 //p(Medium | evidence) 

0.0; //pfLow | evidence); 

} 

Posterior distribution; 

probability ( "Greenj3d"){//1 variable(s) and 5 values 
table 

0.0 If P(Very__high | evidence ) 

0.0 //p(High I evidence) 

1.0 //p(Mediumi evidence) 

O.o ;/p(Low | evidence) 

0.0; // p(Very_low l evidence ); 


I E.«~™ I E.IW.fTL * .! - ; - 1 

Fig 6.8 Probability Prediction Obtained for GBD (Medium 1) 
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P(High | evidence) 
//p(Medium | evidence) 
//p(Low| evidence); 


temp 


rushing_str< 


SIS. iBM 

£<ffl Network 


f • 1 1 




iPHWw I wWiwiiP 


Uombuston^chamberjemp 


GrlTnjemp 


No^ofj'olls Mlx_t#mp Storage.jlfn* 


Fomiing - .pr#**ur# 


P^*;_mix_t*m|® n 8- tem P* rir, 8 


Graphite 


pS c fi 

X, C«B»djporotfty ! R 
oxidation resiftaac* JBl ^*8-' 


chimfcaLpurity 


lag^reslstance 


carbon content 


MitiI M powd«r 


Posterior distribution: 

robability ( "Pitch" ) { m variable(s) and 3 values 
table 

1.0 //p(High | evidence) 

0.0 //p(Medium | evidence) 

0.0; Si p(Low | evidence ); 


osterior distribution: 

robability ( ,, CokedjDorosity"){//1 variable(s) and 3 values 
table 

1.0 // p(high 1 evidence ) 

0.0 //p(medium | evidence) 

0.0; tfp(low| evidence); 


osterior distribution: 

robability { "Coked.porosity'’ ) { //I variable(s) and 3 values 
table 

1-0 //p(high | evidence) 

0.0 // p(medium | evidence ) 

0.0; //p(low| evidence); 

osterior distribution: 

robability ( "Green.bd" ) { m variable(s) and 5 values 
table 

0.0 // p(Very_high I evidence ) 

0.0 //p(High | evidence ) 

1.0 // p(Medium | evidence) 

0.0 // p(Low| evidence) 

0.0; // p(Very_low | evidence ); 





■Jevefteyvt I *blM»§ 


6.9 Probability Prediction Obtained for Coked 


create 


robability ( "Qold_crushing__strength , ’){//l variable(s) and 
table 

0.0 // PCVen,c„high | evidence ) 

1.0 ;/p(High | evidence) 

0.0 ;/p(Medium | evidence ) 

0.0 //p(Low | evidence) 

0.0; ff p(VeryJow | evidence ); 


heater 


■umjtmp 


Retention time 


Combusfon_chamber__temp 


osterior distribution: 

robability ( ,, carbon_contenf ) { //I variable(s) and 5 values 
table 

0.0 Si p(very_bigh | evidence ) 

1.0 //p(high | evidence) 

0.0 //p(medium | evidence) 

0.0 // p(low 1 evidence ) 

0.0; //p(ver/_lowj evidence ); 


Mixmg^time Gn«n_»mp 


No_©fj , ©«* Mix^temp StoFage^tlme 


Jtingjemperirij 


Foqming^pressure Press^mix^tem 


'osterior distribution: 


robability ( "Metal_powder"){//1 variable(s) and 3 values 
table 

0.0 //p(high | evidence) 

1.0 //pimedium | evidence) 

0.0; //p(low | evidence); 


Graphite 


Cofd_crushing_str< 


osterior distribution: 

robability ( "HMOR" ) { in variable(s) and 6 values 
table 

1.0 //p(WH I evidence) 
0.0 //p(VHl evidence) 

0.0 //p(H | evidence) 

0.0 //p(M| evidence) 

0.0 // p(L | evidence ) 

0.0; // p(VL | evidence); 


ColKed^porosity 


:|ag_resistance 


oxTdation_.resi stance 


carbonj:ontent 


powder 


d.t Vaiahte | 6dlt Fi.ir.dio n | Edit Nety^vork — ■' — -7 ~ r T 

Fig 6. 1 0 Probability Prediction Obtained for HMOR (Very-Very- ig ) 
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'mv,'! 

itij * i $N*t t I >h m 


rriESj^voB ayes Console 


Ofosor ve /** °P tlons Help 

I To query on a particular node, click on it. 


mm 


mmm 


Plf^r|pw]l'l<ivoBdyes Console 




. lit 


Crtalt I Mm** I I Qu$ff 


m 

t GrSm fc h«it«rJriimjemp 
itfonjim© 

Combusion^chamberj^mp 


mTngjxm* Grain j#mp 


Mix#rjt#mp 


No w of_roJU Mix^t#mp Storagt^time 


F^ming_pr«s»ur© Pr#ss_mix J;@m^ x *® tt * ns - t@mper ' ,n ^ , 


Graphit© 


Grttti fed 


CoTd^crushing^stn 


C©K©<i.j><>r<»»lty 


oxidation r©ti»t»nc© 


cb©micai_purity 


>lag_r«slstanct 


carbon content 


M©tiJ,j>©wd#r 


Posterior distribution: 

probability ( "carbon_contenf ){//i variable(s) and 5 values 
table 

0.0 // p(very_high | evidence ) 

1-0 //p(high | evidence) 

0.0 //p(medium | evidence) 

0.0 //p(low | evidence) 

j °-0; //p(veryJow | evidence); 

Posterior distribution: 

probability ( "MetaLpowder" ) { //I variable(s) and 3 values 
table 

0.0 //p(high | evidence) 

1.0 //p(medium | evidence) 

j 0.0; //p(low | evidence); 

Posterior distribution: 

probability ( "Green_bd" ) { //I variable(s) and 5 values 
table 

0.0 //p(Very_high [ evidence ) 

0.0 //p(High | evidence ) 

1.0 //p(Medium | evidence) 

0.0 //p(Low| evidence) 

0.0; //p(VeryJow | evidence); 

} 

Posterior distribution: 

probability ( "oxidation_resistance"){//1 variable(s) and 5va 
table 

0.0 ff p(very_high | evidence) 

1.0 //p(high | evidence ) 

0.0 //p(medium | evidence) 

0.0 //p(low[ evidence ) 

0.0; //p(veryJow | evidence); 


• , «! 5 5 S li.S«i*i t 

fiMftJKMA | BdllNolwork 


Fig 6.1 1 Probability Prediction Obtained for Oxidation Resistance (High=l) 


t Griin.hiat^rJrumjimp 

itton tlm# t . . 

Combu$ion_chamberj:emp 


f^tx©r - t«mp 


GraTn_t#mp 


No^of^roliis Storage^tim# 

F J|ing_ P r*.,ur* M - m.x_«m^ n »- t * mp#rint 


GrSphlt© 


jporostty 


CpT3_c 

c^mi 
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Fig 6.12 Probability Prediction Obtained for Slag Resistance (High 1) 
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FIG 7.23 Surface Plot for Mixer Temperature & Mixing Time Vs Mix Temperature 
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FIG 7.30 Plot for Graphite Vs GBD 
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FIG 7.34 Plot for Graphite Vs GBD 
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Appendix I 


RESULTS OBTAINED ON SAMPLE DATA APPLIED TO THE FIS 

SYSTEMS 


Table 7.11 Results for Grain Temperature 
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Table 7.12 Results for Mix Tempera ture 
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Table 7.13 Results for PMT 
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Table 7.13 Results for GBD 
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Table 7.14 Results for COKED POROSITY 
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Table 7.15 Results for CCS 
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Appendix J 


Table 8.4 Results for 86 data set obtained using Matlab (for trimf & trapf) Vs 
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is MEDIUM) and (ROLL is HIGH) and (FP is HIGH) then (GBD is out1mf234) (1) ‘ 
is HIGH] and (ROLL is LOW] and (FP is LOW) then (GBD is out1mf235) (1) 
is HIGH) and (ROLL is LOW) and (FP is MEDIUM) then (GBD is out1mf236) (1) 
js HIGH) and (ROLL is LOW) and (FP is HIGH) then (GBD is outl mf23?) (1) 
is HIGH) and (ROLL is MEDIUM) and (FP is LOW) then (GBD is outl mf238) (1 ) 
is HIGH] and (ROLL is MEDIUM) and (FP is MEDIUM) then (GBD is outl mf233) (1 ) 
is HIGH) and (ROLL is MEDIUM) and (FP is HIGH) then (GBD is out1mf240] (1) 
is HIGH) and (ROLL is HIGH) and (FP is LOW) then (GBD is outl mf 241 ) (1)‘ 
is HIGH) and (ROLL is HIGH) and (FP is MEDIUM) then (GBD is out1mf242) (1) 
s HIGH] and (ROLL is HIGH) and (FP is HIGH] then (GBD is out1mf243) (1) 

and and 


Connection 
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